



THE UNIVERSITY OF CHICAGO 


ESTIMATION OF VENUS WIND VELOCITIES FROM 
HIGH-RESOLUTION INFRARED SPECTRA 


A DISSERTATION SUBMITTED TO 
THE FACULTY OF THE DIVISION OF THE PHYSICAL SCIENCES 
IN CANDIDACY FOR THE DEGREE OF 
DOCTOR OF PHILOSOPHY 

DEPARTMENT OF THE GEOPHYSICAL SCIENCES 


BY 

MftRY ANN HRIVNAK SMITH 


CHICAGO^ ILLINOIS 
DECEMBER 1978 



acknowledgements 


The author would like to e:^ress her sincere thanks to Professor 
Lewis D- Kaplan for his guidance, encotiragement and support for the work 
leading to this dissertation. 

Special words of appreciation are also necessary for Dr. Pierre 
Connes and Dr. Janine Connes, of the Centre National de la Eecherche 
Scientifique, Orsay, who measiired and processed the planetary spectra and 
sv 5 »plied the con^uter program and other data necessary for the detailed 
stu<fy of these spectra. Thanks are also due to Dr. Guy Guelachvxli of 
the Laboratoire d’infrarouge, university of Paris, Orsay, who on request 
made laboratory measurements of the pressure-induced lineshifts of HCi and 
HF. Without these data this 'dissertation could not exist. 

The author also thanks C. Boulet of the Laboratoire d’infrarouge 
for providing his calculations of the lineshifts for ECU perturbed by HF 
and for his helpful comments on the work disc\issed in Chapter II. She is 
also grateful to Professor Gilbert Amat of the Laboratoire de Physique 
Moleculaire et d'Optique Atinospherique , C.N.R.S., Orsay, for his helpful 
discussions concerning this work, and to Professor Dave Fultz of the Uni- 
versity of Chicago for his comments on the results presented in Chapter IV. 

Acknowledgement is also due to Mrs. Mary Adlyn Evans for her typ- 
ing of this dissertation.' 

This research was supported by the U.S. National Aeronautics and 
Space Administration, under grant NGR 14-001-194. 

Finally, the author would like to e:q)ress her deep appreciation to 

ii 



her hiasband, Charles, and to her mother, brothers and her late father for 
their past and continued support and encoiaragement. 


iii 



TABIiE OF CONTENTS 


ACKNOWLEDGEMENTS ii 

LIST OF TABLES vi 

LIST OF ILLUSTRATIONS vii 

ABSTRACT X 

Chapter 

I. INTRODUCTION 1 

II. PRESSURE-INDUCED SHIFTS OF INFRARED ABSORPTION LINES 3 

2.1 Characteristics of the shift 

2.2 Theoretical discussion 

2.3 Pressure- induced lineshifts for molecules 
common in planetary atmospheres 

III. THE VENUS ATMOSPHERE 13 

3.1 Composition and structure 

3 . 2 Atmospheric motions 

IV. ESTIMATION OP VENUS ZONAL VELOCITY PROFILES USING 

PRESSURE-INDUCED LINESHIFTS 28 

4.1 Despcription of the spectra 

4.2 Geometry of the problem 

4.2.1 Effect of translational and rotational 
Doppler shifts 

4.2.2 Scaling of residual Doppler shifts to 
obtain Venus equatorial zonal winds 

4.3 The direct problem: calculation of expected 

line asymmetries for a model atmosphere 

4.4 The inverse problem: estimation of the velocity 

profile from measured line asymmetries in the HCJl 
2-0 band 


4.4.1 Procedure 

4.4.2 Results and error analysis 



TABLE OF CONTENTS — Continued 


Chapter 

V. CONCLUSION 63 

EEFEBENCES 113 


V 



LIST OF TABLES 


Table Page 

1. Physical characteristics of Venus and Earth 66 

2. Composition of the Venus atmosphere 67 

3. Parameters for recorded spectra of Venus 71 

4. Hydrogen halide lines observed in the Venus spectra ..... 72 

5. Venus Doppler shifts calculated from Ephemeras 75 

6. Reference heights of chords across spectral 

lines in the A-H(l) model 76 

7. Venus wind velocities calculated from CO 2-0 band 77 



LIST OF ILLUSTEATIOMS 


Figure Page 

1. Measured pressure-induced lineshifts for HP per atm. 

Ar. (a) 1-0 jaand. Circles: Guelachvili and Smith 

C1978) ; Triangles; Jaffe ^ al. (1965) ; Crosses: 

Oksengorn (1963). (b) 2-0 band. Circles: Guelachvili 

and Smith (1978); Crosses: Wiggins et al. (1970) 78 

2. Measured lineshifts per atm at in the HCil 2-0 band. 

(a) H^®CA. Circles: Guelachvili and Smith (1978) ; 

Triangles: Eank et al. (1960); Inverted triangles: 

Rank et al. (1964); Crosses: Hirshfeld et al. (1960) 

(b) H^Cil. Circles; Guelachvili and Smith (1978); 

Crosses; Hirshfeld et al. (I960) 79 

3. Measured lineshifts per atm CO 2 in the ECS, 2-0 band. 

(Guelachvili and Smith/ 1978). (a) (b) H^^Cl ... 80 

4. Self-induced lineshifts in the HC£ 2-0 band measured 

at low pressure, circles: Triangles: h^^c£. 

(Guelachvili and Smith/ 1978) . Dashed line: Kimel 

et al . (1959) 81 

5. Self-induced lineshifts in the HP 2-0~band measured at 
low pressure. Circles; Observed shifts (Guelachvili 
and Smith, 1978). Dashed line; Shifts calculated by 

Boulet et al. (1976) 82 

6. Ground-state rotational energy levels and near-resonant 

transitions for H^^Ci and HP (after Boiilet et al. 1976) . . 83 

7. Lineshifts in the HCX. 2-o band perturbed by 10 torr HP. 

Circles: H^®c£ observed shifts; triangles: observed 

shifts (Guelachvili and Smith, 1978); inverted triangles: 
calculated shifts (Boulet, 1977) 84 

8. Lineshifts in the HP 2-0 band perturbed by 10 torr HC5- 

(Guelachvili and smith, 1978) 85 

9. Measured lineshifts per atm n 2 in the HCJl 2-0 band. 

Circles: Guelachvili and Smith (1978); triangles: Rank 

^ . (1960a); Crosses: Shifts measured by Kimel et al. 

(1959) for 1 atm dry air. (a) H^®C&. (b) H^^C£ 86 


vii 



LIST OP ILLUSTRATIONS-- Contxnued 

Figure Page 

10. Measured lineshifts for HP perturbed by COa. (a) 

1-0 band. Circles; Guelachvili and Smith (1978) ; 

Triangles; Shaw and Lovell (1969) , (b) 2-0 band. 

(Guelachivili and Smith, 1978) 87 

11- Profile of temperatures in the Venus atmosphere 88 

12. Measurements of Venus zonal winds 89 

13. Distribution with latitude of Venus zonal veloc- 
ities - 90 

14. Interferometric spectra of Venus near 5750 cm 

(a) Venus A. (b) Venus 5 91 

15. Residual lineshifts from Venus spectra for the 

HCil 2-0 band., (x lo"^ cm"^) (a) (b) 

93 

16. Residual lineshifts for HP in the Venus spectra ' 

(x 10 ^ cm ^). (a) 1-0 band. (b) 2-0 band 94 

17- Schematic drawing of calculation of the asymmetries 

in spectral lines by the DERPTH program 95 

18. Vciriation of translational Doppler shift factor from 

1971 to 1980 (Niehaus and Petrie, 1961) 96 

19. Appearance of the Venus disk at the time of observa- 
tions. (a) Venus A. (b) Venus 5 97 

20. Variatxon of brightness along parallels on the planetary 
disk. Solid curves; Lambert's Law reflection. Dashed 

Guinres; isotropic scattering 98 

21. Relative contribution to the earthward component of 

Venus zonal velocity. (a) Venus A. (b) Venus 5 99 

22. Pressure, temperature and velocity profiles for Venus 

model atmospheres 100 

23. Absorption coefficient for significant levels in the 

median model atmosphere for the R4 line. (a) 

Venus A. (b) Venus 5 101 


viii 



LIST OF _ ILLUSTRATIONS — Continued 


Figure Page 

24 . Transmittance between 95 km and significant levels 

in the median model atmosphere for the R4 

line. (a) Venus A- (b) Venus 5 103 

25. Calculated asymmetric profile for the R4 line 

in the medxan model atmosphere for Venus 5 105 

26. Calculated lineshifts in the H^^CA R-branch for line 
profile "D” in the A-H{1) model atmosphere. (a) Venus 

A. (b) Venus 5 106 

27. Calculated lineshifts in the H^^cA R-branch for line 
profile "D" in the median model atmosphere, (a) Venus A. 

(b) Venus 5 107 

28. Calculated lineshifts in the H^^CA R-branch for line 
profile "D" in the A-H(2) model atmosphere. (a) Venus A. 

(b) Venus 5 108 

29. V Zonal velocity profiles retrieved from model results. 

Triangles: Venus A; Circles: Venus 5; Dashed lines; 

Original velocity profiles 109 


30. Weighted-average lineshifts in the HCA 2-0 band for the 

averaged planetary spectra. (a) Venus A. (b) Venus 5 . . . 110 

31. Venus zonal velocity profiles estimated from inter- 
ferometric spectra. Triangles: Venus A; Circles; 

Venus 143 (first day); Rhombi: Venus 161 (last day). 

Diamonds: Venus 5; Squares; Venus 226 (first day) ; 

Crosses: Venus 234 (last day) ..... Ill 

32. Longitudinal cross section of Venus equatorial zonal 

velocities. Dashed lines: isotachs; contour interval 

is 50 m sec ^ . Solid lines ; estimated streamfunction 

in units of m sec ^ bar 112 


ix 



ABSTRACT 


--Zonal velocity profiles in the Venus atmosphere above the clouds 
Cein be estimated from measured asymmetries of HCJl and HF infrared absorp- 
tion lines in high-resolution Fourier interferometer spectra of the planet. 
These asymmetries are caused by both pressure-induced shifts in the posi- 
tions of the hydrogen-halide lines perturbed by CO 2 and Doppler shifts due 
to_ atmospheric motions. Particularly in the case of the HCA 2-0 band, the 
effects of the two types of lineshifts can be easily isolated, making it 
possible to estimate ■ a profile of average Venus equatorial zonal velocxty 
as a function of pressure in the region roughly 60 to 70 km above the sur- 
face of the planet. Although there are large uncertainties in the result- 
ing wind profiles, the mean magnitudes of the winds appear to be realistic 
and are not inconsistent with those measured by other investigators. The 
mean profiles obtained’ show strong vertical shear in the Venus zonal winds 
near the cloud- top level, and both the magnitude and direction of winds 
at all levels in this region appear to vary greatly with longitude relative 
to the sub- solar point. 
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CHAFTER I 


INTRODUCTION 

The development of a high-resolution Fourier-transfona infrared 
interferometer in 1966 by J. and P. Connes has been a significant contri- 
bution to the study of planetary atmospheres. The resolving power of the 
Fourier-transfozm interferometer in the infrared spectrum of Venus was 
about one hundred times greater than that of the, best scanning grating 
spectrometers in use at the time (Connes and Connes , 1966) . 'A major advan- 
tage of Fourier spectroscopy is that line positions can be measured with 
great accuracy, so that positions determined in planetary spectra often 
approach the accuracy of laboratoiry measurements. Details of the design 
of the Connes instrument and the Fourier-transform scheme used to produce 
the final spectrum can be found in several of their papers (Connes and 
Connes, 1966; Connes et al., 1967). 

The earlier Connes planetary spectra, published in their Atlas in 
1969, had a resolution of 0.08 cm which was far superior to that of pre- 
vious planetary spectra. Line positions in the HCi 2-0 band in the Venus 
spectra were measured with an accuracy better than 0.01 cm Prom these 
planetary spectra the Connes and their associates were able to identify 
and measure for the first time the presence of CO, HCil and HF in the Venus 
atmosphere (Connes ^ , 1967, 1968) and CO in the Martian atmosphere 

(Kaplan et a^. , 1969) . 

During October and November of 1973, using an improved version of 
the interferometer. Dr. P. Connes took interferograms of Venus and Mars 
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through the Hale 5in‘ (200in) telescope at Mt. Palomar. A second set of in- 
terferograms of Venus were taken at the same site in June 1974. In addi- 
tion a lunar spectrum was measured during the 1973 observation period, for 
comparison and to aid in identifying telluric absorption features in the 
planetary spectra. The resolution of all the spectra is about 0.01 cm 
and line positions can be determined to the necirest 0.001 cm ^ in the 1.7 
ym region. The Mars spectra measured in 1973 have been used to provide 
a pre-viking estimate of the amount of inert gases present in the Martian 
atmosphere (Kaplan et al . , 1975) . 

The accuracy of line positions measured in these new planetary 
spectra is sufficient to allow detection of pressure-induced shifts of HCit 
and HF infrared absorption lines due to the C0_ in the Venus atmosphere. 

A 

In addition to shifts measured at the peaks of the absorption lines, there 
were found asymmetries in the line shapes apparently arising from both the 
pressure-induced lineshifts and Doppler shifts due to Venus atmospheric 
motions. It is possible to invert these measured line asymmetries for a 
given absorption band to obtain an estimated profile of the average equa- 
torial zonal velocity as a function of pressure for the Venus atmosphere 
above the main cloud deck. 

Before the estimation of Venus wind profiles from absorption line 
asymmetries is described, it is very helpful to discuss in more detail the 
phenomenon of pressure-induced lineshifts. 



CHliPTER II 


PRESSOEE-INDUCED SHIFTS OF INFEaRED ABSORPTION LINES 

A detailed and precise understanding of the vibration-rotation 
bands of the hydrogen halides is important for high-resolution infrared 
spectroscopic studies of the atmospheres of both Earth and Venus. Since 
these halogen compounds, primarily HCZ and HF, are at most minor consti- 
tuents of the planetary atmospheres, their absorption spectra are very 
likely to be affected by other atmospheric gases. The most common effect 
is Lorentz, or pressure- induced, broadening of the absorption line. How- 
ever, in addition to the broadening of the line profile, the line center 
is displaced by a small amount proportional to the pressure of the per- 
turbing gas. 

Laboratory measurements of pressure-induced lineshifts for hydro- 
gen halides were first made in the late 1950s (Kuipers, 1958; Kimel et al. , 
1959; Rank et al. , 1960a, b) . The original intent of these early experi- 
ments was to measure foreign-gas broadening and self-broadening in the HCA 
and HF 1-0 and 2-0 bands, using gas pressures from several cmHg to 5 atm. 
Specific measurements of pressure-induced lineshifts, as well as theoreti- 
cal studies, have continued to the present time. The most commonly studied 
lineshifts have been those of HF or HC5- perturbed by the noble gases 
{Boulet ^ al . , 1972a, b; Jarecki and Herman, 1975) . More recently theo- 
retical studies have appeared dealing with pressure-induced lineshifts due 
to more polar perturbers (Boulet et al., 1976). 
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2.1 Characteristics of the Shifts 

The most important characteristic of pressure-induced lineshifts 
in terms of their usefulness in the study of planetary atmospheres is that 
the magnitude of the shift for a given spectral line varies linearly .with 
the pressure of the perturbing gas. This linear relationship became obvi- 
ous in experiments by Kimel ^ a^. (1959) , where self-induced lineshifts 
in the HCil 2-0 band were measured at several gas pressures between 20 and 
700 torr, and the linearity was verified for other molecules by later in- 
vestigators (Rank et al. / 1960a; Shaw and Lovell, 1968), If a set of pres- 
sure-induced lineshifts for a specific pair of absorbing and perturbing 
gases is measured in a planetaiy spectrum, the atmospheric pressure at which 
the maximum absorption occurs could be easily determined from the linear 
dependence of the lineshifts on pressure. 

When the pressure-induced lineshifts in the same vibration-rotation 
band are examined together, it becomes apparent that the lineshifts are 
strongly dependent on the rotational quantum number (J) of the absorbing 
molecule. For example, in Figure 1(a) for the HP 1-0 band perturbed by 
flr, the shifts become increasingly negative for larger values of J, reach- 
ing a maximum magnitude for J values of 6 or 7, then slightly decreasing. 
Pressure— induced lineshifts also show some dependence on the vibrational 
quantum number of the absorbing molecule. This fact can also be seen in 
Figure 1 (b) , where the lineshifts for the HF 2-0 band perturbed by Ar ap- 
pear to be approximately twice as large as corresponding lineshifts in the 
HF 1-0 band. 

The J-dependent pattern of the lineshifts is specifically tied to 
the particular absorber-perturber pair. Rank et al. (1960b) have shown 
empirically for HCA absorption lines that the magnitude of the lineshift 
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depends largely on the polarizahility and molecular mass of the perturbing 
gas. In addition it has been observed that nonpolar molecules such as 
and the noble gases produce smooth, father symmetrical patterns when per- 
turbing strongly polar molecules such as HF and HCX. (Figures 1 and 2) , 
However, the pattern of lineshifts becomes increasingly complex for more 
polar perturbing gases (Figures 3-5, 7, 8) due to stronger electrostatic 
interactions between the absorbing and perturbing molecules. 


2.2 Theoretical Discussion 

' - Broadening effects in the absorption spectra of diatomic polar 

molecules have long been a subject of theoretical investigation, but only 
recently has similar attention been focused on the corresponding pressure- 
induced lineshifts. Most theoretical studies of lineshifts have dealt 
with interactions of diatomic polar molecules with the noble gases (Herman 
and^Tipping, 1970; Boulet et al ■ , 1973a, b; Jarecki and Herman, 1975). Very 
recently Boulet et al . (1976) have published a theoretical analysis of 
pressure-induced lineshifts due to the interaction of several species of 
diatomic polar molecules with each other. The theory of pressure-induced 
lineshifts will be qualitatively summarized here, following the model and 
notation used by the several authors mentioned above. 

The foreign-gas- induced shift of any vibration-rotation line may 
be considered as the difference of the shifts of the two energy levels in- 
volved in the transition, that is. 


(S 


line 




( 2 . 1 ) 


where the subscripts i and f denote the initial and final levels, respec- 
tively. Each level shift has both a first-order and a second-order con- 


tribution 
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where and are the vibrational and rotational quantum numbers for 
the level in the absorbing molecule, and v^ and J^~are the corresponding 
quantum numbers describing the rotation-vibration state of the perturbing 
molecule. Of cotirse v^ and are not applicable to monatomic perturbers. 
The two terms in the level shift arise from the expansion of the 
binary intermo'lecular potential V using the impact approxxmation developed 
by Foley (1946) and later applied specifically to the theory of pressure- 
induced lineshifts by Herman and Tipping (1970), Boulet et al. (1972a,b) ; 
Jarecki and Herman (1975) and Boulet et al. (1976) . The first-order 


term in the expansion of V is imaginary and contributes only to shifts of 
spectral lines. There are a number of second-order terms which are com- 
plex; the real part of each of these terms- contributes to line broadening 
and the imaginary part to the lineshifts. 

In the calculation of lineshifts, each term in the expansion of 
the binary intemolecular potential must be averaged over all possible 
directions of the initial relative velocity of the colliding molecules in 
order to obtain the net perturbation of the level in the absorbing mole- 
cule. As a result of this averaging, the first-order term vanishes when 
the absorbing molecule is in the ground vibrational state, i-e. , = 0 

for v^ = 0. For higher vibrational states, the first-order level shift 
is always negative and less important than the second-order shift except 
for very large values of v^ and J^. The magnitude of the first-order level 
shift increases directly with v^, so that always pro- 


ducing a negative (or zero for purely rotational lines) first-order 



7 


contribution to the line shift. 

The- imaginary parts of the second-order terms in the collision 
Hamiltonian (derived from the e:^ansion of the binary intermolecuiar poten- 
tial V) provide the most significant contribution to the pressure-induced 

I 

lineshifts. The predominant shifting mechanisms vary according to the im- 
pact parameter b, which represents the effective collision cross-section 
of the absorbing molecule relative to the perturber. Below a certain 
critical "cutoff value" b^, short-range interactions predominate, as is 
the case for diatomic polar molecules perturbed by noble gases (Herman and 
Tipping, 1970; Jarecki and Herman, 1975). Among the important short-range 
shifting mechanisms are anisotropic repulsion forces which contribute 
mainly to the shifts for low J values by causing bending of the particle 
tra.jectories. Both elastic and inelastic collisions produce phase shifts 
in the vibration-rotation mode of the absorbing molecule. Long-range elec- 
trostatic interactions may be safely neglected. 

For values of the impact parameter much greater than b , as is the 

"o 

case for diatomic polar molecules as perturbers, the interactions between 
the absorbing and perturbing molecules are primarily long-range and elec- 
trostatic. The variety of short-range forces which were important for noble- 
gas perturbers may be ignored in this case. Long-range dipole-dipole in- 
teractions cause the colliding molecules to induce rotational transitions 
in each other. These induced transitions occur for all combinations of 
absorbing and perturbing polar molecules, but in cases where the rotational 
constant for the perturber is approximately twice that for the absorber (as 
for HCA-HP, DC5.-HC£ and DF-HF) , near-resonant induced transitions occur 
which are particularly efficient in causing level shifts. As illustrated 
in Figure 6 for HC5r perturbed by HF, there are two types of these near- 
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resonant interact ions - 

(i) In the first type (I) of near-resonant interaction, the transi- 
tions in resonance differ in energy by approximately integral multiples 

of the ground- state, rotational constant for HF. Since for HF the 

OHF 

-7 

upper vibration-rotation states have populations that are at most 10 
the population of the ground state at normal atmospheric temperatures, 
only the ground-state rotational transitions in HF may be considered as 
effectively producing level shifts in HCS, molecules. Several of these 
Type I resonances for both even and odd levels are denoted by the dashed 
arrows, e.g., Al< — >a2, Bl< — >B2, m Figure 6. The most effective of these 
resonances in producing HC£ level shifts are those whose energy difference 
Aoi is least. For a first approximation contributions from resonances with 
energy differences larger than ± can be neglected. 

The largest part of the Type I contribution to the,HC£ level shift 
is represented by the difference of two terms corresponding to the Al< — >A2 
resonance and the Bl< — >B2 resonance. Each term contains a Boltzmann 
factor expressing the population of the initial level of the perturbing 
HF transition: exp[{-hcAT) *B (J'.+l) }, where J'. = f(J.) is a rela- 

tionship defined by the two transitions in resonance. In each term this 
Boltzmann population factor is also multiplied by an "interruption function" 
IFl(Au), where the energy difference Aw = B ± Ab (J.+l) for even J., 

OHF VI i 

or AoJ = B.„_ ± Ab Cr. for odd J. , with AB = 2B - B The function 

OHF VI 1 V vHCJl OHF 

IFl describes the effectiveness of a given pair of transitions in producing 
level shifts in the absorbing molecule. For HCA and HF at typical atmo- 
spheric temperatures, IFl (Aw) has a maximum at Aw = 5.7 cm ^ and is linear 
for smaller values of Aw and inversely proportional for larger values. 

Thus the sign and magnitude of the Type I contribution to the HCil level 
shift depend on the signs and magnitudes of the Ab^ and the magnitude of 
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B [through IFl(Au)], as well as on the relative populations of the HF 
OHF 

states involved in the two primary pairs of near-resonant transitions. 

(ii) The second type (II) of shift-producing resonance involves 
transitions which have approximately the same energy, i.e., Au = for 

even J. and Aw = AB (J.+l) for odd J.. ‘The Type II resonances are de- 
noted by small letters, e.g., a, b, c, in Figure 6. As is obvious from 
the figure, in the ground state of HCi the resonances cause the odd-J 
levels to be perturbed downward and the even-J levels to be perturbed up- 
ward. In the same manner as for the Type I resonances, the Type II con- 
tribution to the level shift is proportional to the population of the 
lower HP state involved xn the resonant, transition and to the interruption 
function IFl{Aw) , although for the Type II level shift there is only one 
pair of resonant transitions being considered, since in thxs case Aw is 
directly proportional to Ab^, the sigh of Ab^ is particularly important. 

The level shifts will be positive for even and negative for odd if 
Ab^ is positive (as it is for the ground state of HCi) , but if Ab^ is nega- 
tive (as in the second vibrational state of HC£) the phase of the shift 
oscillations will be reversed. 

As illustrated by Boulet ^ al . (1976) , the magnitudes of Ab^ and 
B^ for the perturber determine the relative influence of Type I and Type 
II contributions to the shifts for a given vibrational state of the absorb- 
ing molecule. In the case of DCil perturbed by HCS,, the relative influence 
of Types I and II resonances is approximately the same in the first and 
second vibrational states, while in the ground state the Type II resonances 
clearly predominate. Moreover, the magnitudes of the ground-state level 
shifts are generally larger than those in the upper states. When line shifts 
for the DCA 2-0 band are calculated, the predominance of the- Type II 
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ground-state shifts results in a strongly oscillating pattern of line 
shifts, corresponding to those observed in the Iciboratory by Jaffe, et al. 
(1964)', 

For HC& perturbed by HF, the values for all the lower vibra- 
tionah states are large enough relative to so that the Type II reso- 

nances are the predoininant shifting mechanism in all these states. The 

magnitude of AB 2 is about three times as large as that of Ab , making the 

o 

shifts in the second vibrational state generally larger than those in the 

ground state; however, since Ab is positive and ^02 negative, the shift 

o 

oscillations in the two states are out of phase. The resulting' line shifts 

1 

for the 2-0 band of HCA, both calculated (Boulet, 1977}* and observed 
(Guelachvili and Smith, 1978} , are shown in Figure 7. 

The near-resonant interactions of ucAand HF can also be used to 
predict lineshifts for HF vibration-rotation transitions perturbed by HC£. 
The calculation of lineshifts becomes more complex sxnce each HF rotational 
transition is involved xn four major Type I near-resonant interactions with 
ground-state HCfi. rotational transitions and in two opposing (but not equal} 
Type il resonances (Figure 6}. In addition, both isotopic species 
and are interacting with HF molecules in proportion to their rela- 

tive abundance. Thus the number of second-order terms in this calculation 
is quadruple the number necessary for the calculation of HCZ shifts due 
to HF. A partial calculation of the net relative second-order contribu- 
tion to the 'HF lineshifts by HCiJ. seems to correctly predict the behavior 
for low J values of the HCfi.-induced lineshifts observed in the HP 2-0 band 
by Guelachvili and Smith {1978} , shown in Figure 8. However, more detailed 
theoretical calculations of these lineshifts are necessary. Nevertheless, 
it can be observed from Figures 7 ‘and 8 that shifts of RCZ lines perturbed 
by HF possess a strongly oscillating character, while for HF perturbed by 
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HCS. there is no real obvious oscillation. 

The observed pressure-induced lineshifts of ECS, eind HF shown in 
Figiares 4, 5, 7 and 8 were deduced from experiments which were not performed 
for that purpose. The gas pressures used (not more than 10 torr) are not 
high enough to get very accurate measurements of the shifts. The results 
are' nevertheless .significant, since they contain the only HC5.-HF line- 
shift measvtrements in existence at the present time. Obviously more theo- 
retical and experimental study of this problem is necessary. 

2.3 Pressure-induced Lineshifts for 
Infrared-absorbing Molecules in Plane- 
tary Atmospheres 

For application to the study of , planetary atmospheres, we are pri- 
marily interested in pressure-induced lineshifts of trace-constituent in- 
frared-absorbing gases perturbed by the major atmospheric constituents. 

For the earth's atmosphere, measurements of CO, CO 2 / H 2 O, O 3 , N 2 O, CH 4 , SO;., 

NO, NO 2 , NH 3 , HC£, HF or HNO 3 perturbed by N 2 , O 2 and Ar would all be im- 
portant. Unfortunately, none of these shifts have been measured in the 
laboratory with the exception of CO perturbed by Ar (Ranic et al . , 1963 ; 
Bouanich, 1972) and HC^. and HF perturbed by N 2 and by Ar (Rank ^ al . , '1960a; 
Rank et al. , 1963; Hirshfeld et al . , 1960; Jaffa et al. , 1965; Oksengorn, 

1963; Wiggins et al . , 1970; Guelachvili and Smith, 1978) . In Figure 9 
measurements by Guelachvili and Smith (1978) for the HC£ 2-0 band perturbed 
by N 2 are compared with measurements by Kimel ^ ^1. (1959) for the same 

band perturbed by 1 atm dry air. Since dry air contains 78 percent N 2 by 
volume, it is not unexpected that the two sets of lineshifts are quite simi- 
lar. Note also that the nitrogen curves are similar in form and magnitude 

I 

to the ones for argon (Figure 2) ; this is not surprising since both mole- 
cules are nonpolar and have roughly similar molecular weights. The differences 
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beiween the shifts by air and by nitrogen can probably be accounted for 
by the effects of the 21 percent O 2 and 0.9 percent Ar also present in 
dry=^air. unfortunately, there exists only a single measurement of a line- 
shli^ft-in the HCi 2-0 band perturbed by O 2 ; Rank et al. (1960a) have measured 
that the RIO line of is shifted -22.0 x 10 ^ cm ^ per atm O 2 . With 

sb^little information available, meaningful quantitative compcirisons for 
the- earth's atmosphere cannot yet be made. 

por the Venus atmosphere the primary interest is measured shifts 
of HCA and HF perturbed by CO 2 , although measurements of CO perturbed by 
CO? -wo uld be helpful if they existed. It is interesting to note that the 
pressure- induced shift patterns for HC£ and HP perturbed by CO 2 (Figures 
3 and 10) contain more "kinks" and are distinctly different frcrni the cor- 
responding shift patterns for these bands perturbed by Ar (Figures 1 and 
2) . In terms of the theoretical model discussed in the preceding section, 
the shifts by Ar represent mainly the weak first-order contribution plus 
stronger second-order contributions arising from short-range anisotropic 
repulsion forces. For the shifts by C02f the first-order and short-range 
second-order contributions remain, but there are additional long-range 
second-order terms arising from dipole-dipole interactions between the 
hydrogen halide molecule and CO 2 molecules in asymmetric vibrational states. 
Since CO 2 is at best a temporary dipole, the second-order long-range inter- 
actions are not nearly as strong as those between HF and HCi, and strongly 
oscillating shift patterns are not expected to occur. The interactions 
between HCJl or HF and CO 2 have not yet been theoretically studied in de- 
tail, probably because of the complexity of calculations necessary in order 
to include both short- and long-range second-order effects. 



CHAPTER III 


THE VENUS ATMOSPHERE 

Although scxentific interest in the planet Venus began with tele- 
scopic observations by Galileo in 1609, the determination that Venus pos- 
sesses an atmosphere was not made until nearly 200 years later. In the 
latter part of the eighteenth century several observers noted that the 
horns of the crescent of Venus extended far’ beyond the semicircle, which 
is not the case for an airless body such as the moon. There also was noted 
a considerable time difference between the predicted and actual occur- 
rences of dichotomy, the point when exactly half the planetary disk is il- 
luminated; this time difference could only be due to the presence of an 
atmosphere. 

3.1 Composition and Structure 

Information about the composition of the Venus atmosphere was not 
obtained until the early twentieth century, when spectroscopic observations 
became possible- At that time Venus and Earth were considered to be iden- 
tically twins , and the Venus atmosphere was expected to- be quite similar to 
the terrestrial one. The early Venus spectra revealed that the major op- 
tically active constituent of that planet’s atmosphere is ^carbon dioxide, 
not oxygen as had been e:5^ected for a terrestrial-type atmosphere. 

Until the recent successful missions carried out by the Soviet 
probes in the Venera series, studies of the composition and structure of 
the Venus atmosphere were limited to the region above the optically thick, 
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•unbroken cloud deck. Nevertheless, more recent spectroscopic measure- 
ments made at much higher resolution than those at the turn of the century 
have fairly well established the composition of the Venus atmosphere. The 
atmosphere is composed almost entirely of CO 2 , with CO appearing as a uni- 
formly mixed i^nor constituent. The Venus atmosphere also contains small 
but spectroscopically measurable amounts of HCf. and HF (Connes et al., 

1967) , and these two hydrogen halides also appear to be uniformly mixed 
in the atmosphere above the clouds. Despite difficulties due to strong 
tell'uric water vapor absorption features, the amount of H 2 O in the Venus 
upper atmosphere has also been spectroscopically measured. The most re- 
cent of these measurements estimate the amount of water vapor to be roughly 
equivalent to 1 pm precipitable water in the atmosphere above the clouds 
(Gull et a^. , 1974) . Only the five constituents mentioned above have ac- 
tually been detected from spectra of Venus. However, several investiga- 
tors have also calculated upper limits for many other geochemically plau- 
sible minor constituents such as NH 3 / H 2 S, COS, SO 2 and O 3 , based on the 
non-appearance of their absorption features in Venus spectra. - One of the 
nrast comprehensive works of this type is that of Owen and Sagan (1972) , 
who made use of ultraviolet spectra from the Orbiting Astronomical Obser- 
vatory. L. Young (1974) has compiled a table summarizing observations and 
calculated upper limits for the various molecules in the Venus atmosphere; 
the table is reproduced here as Table 2 (a) and (b) , with several more re- 
cent observations added. 

Even though infrared, visible and ultraviolet spectroscopic studies 
were reasonably accurate in estimating amounts of a-fcmospheric constituents, 
estimates of the surface pressure and temperature on Venus varied over a 
very broad range. In the late 1960s and early 1960s radio and microwave 
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astronomy experiments showed that Venus had an unexpectedly high radio- 
brightness temperature in excess of 400°C, Radar experiments also allowed 
investigators to at last determine the true rotation rate of the planet. 
Venus rotates in a retrograde direction, thus differing radically from most 
other planets in the solar system, with a period of about 243 days. Since 
- the orbital period of Venus is somewhat shorter (225 days) , the net. length 
of a "day" from one sunrise to the next is about 117 Earth days. Also, 
the axis of rotation of Venus is nearly perpendicular to the planet's or- 
bital plane, causing a virtual absence of seasons.-. 

Although earth-based microwave measurements had narrowed down the 
wide range of previously estimated temperature and pressure profiles for 
the Venus atmosphere, it remained for the Soviet probes in the Ven^a series 
to provide accurate measurements of temperatures and pressures helow the 
cloud deck. The earlier landers, Veneras 4 through 8, established the 

ki. 

following facts; the surface pressure and temperature are 90-95 atm and . 
740-750 K; the principal atmospheric constituent is CO 2 ; the temperature 
lapse rate in the troposphere (0-50 km) corresponds to adiabatic equilib- 
rium. Veneras 9 and 10 verified the results from the earlier landers. 

In addition, since Veneras 4 through 8 had landed near the morning termina- 
tor and Veneras 9 and 10 landed near the subsolar point, comparison of the 
two sets of ten^erat\ire measurements reveals that daytime temperatures are 
20 to 25°C higher than those at night above the 40-50 km level. 

Direct measurements of Venus atmospheric temperatures and pres- 
sures were made by the Veneras only between the surface and maximum heights 
of '50 to 63 km. Above this region it is necessary to rely on indirect mea- 
surements such as those from the Mariner 5 and 10 flyby missions and from 
earth-based spectroscopic experiments. Figure 11 presents a temperature 
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profile for the Venus atmosphere which is a composite of several of the 
most recent measurements (Marov and Eyabov, 1972; Ainsworth and Herman/ 

1975; Taylor, 1975; Dierenfeldt et al. , 1977; Keldysh, 1977). 

When Venus is viewed from Earth or from space, the visible "sur- 
face" of the planet is formed by an optically thick, unbroken cloud cover. 

In recent years the structiire and composition of these clouds has been the 
subject of considerable study and speculation by planetary scientists. 

However, much of the observational evidence relating to the cloud struc- 
ture and composition is uncertain and sometimes contradictory, and no 
single model consistent with all the observations has yet appeared. 

Significant information about the Venus clouds was obtained by 
Hansen and Hovenier (1974) in their careful analysis of infrared polari- 
zation measurements of reflected light from the planet made by at least 
eight different investigators between 1929 and 1971, From these measure- 
ments they deduced that the Venus cloud particles near the top of the in- 
frared-reflecting layer are spherical in shape (probcibly -liquid droplets) , 
having a remarkably narrow size distribution about the mean particle radius 
of 1.1 pm. Hansen and Hovenier also used these scattering measurements to 
estimate that the cloud-top atmospheric pressure is about 50 mb, which cor- 
responds to a height of about 7,0 km above the surface- This figure corre- 
lates well with pressure of line formation calculated from infrared absorp- 
tion spectra for HC^ (Connes et , 1967) and CO 2 (L. Young, 1972; Dierenfeldt 
, 1977). 

Photographs of Venus made in the visible and infrared portions of 
the spectrum show a rather uniform, yellowish-white cloud "surface." How- 
ever, ultraviolet pictures taken through earth-based telescopes and' by space 
probes such as Mariner 10 clearly reveal light and dark cloud features. 
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According to Keldysh (1977) , these light and dark features correspond with 
high and low areas at the diffuse upper boundary of the clouds. The most 
notable of. the ultraviolet cloud features xnclude a large quasi-permanent 
Y-shaped dark marking which appears to travel around the planet's equator 
with a four-day rotation period, smaller rapidly-moving cloud formations, 
and brightening and darkening of the polar regions (Dollfus, 1975; Limaye,, 

1977) . Since the apparent height of the ultraviolet features was generally 
greater than that estimated for the infrared- scattering cloud, it was specu- 
lated that there might be one or more UV-scattering haze layers above the 
main cloud deck. 

The Soviet probe Venera 8 provided the information that the lower 
boundary of the cloudy region was about 35 km above the planetary surface, 
but gave few clues about the structure or composition of the clouds. A very 
surprising result from this mission was that a considerable amount of visible- 
wavelength sunlight reached the surface of Venus. This encouraged Soviet 
scientists to include television cameras in the instrument packages on 
Veneras 9 and 10, from which were transmitted the first photographs taken 
at the surface of Venus (Keldysh, 1977) . These last two Veneras also car- 
ried a large series of experiments to investigate the optical characteris- 
tics of the atmosphere and clouds. These experiments revealed an apparent 
multilayered cloud structure. The main cloud deck, consisting of a layer 
of particles with radii from 0.5 to 2,5 ym, lies between 49 and 68 km. Be- 
low 49 km there are a few very tenuous haze layers, with no scattering aero- 
sols detectable below 18 km. The upper boundary of the main cloud layer is 
rather diffuse, with considerable variations in particle concentration. As 
mentioned earlier, the light and dark cloud featxires appearing in UV photo- 
graphs seem to be associated with high and low areas in this diffuse cloud top. 
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Although since the Venera 9 and 10 missions the structure of the 
clouds is more certain, their cos®30sition remains a subject of some debate - 
Literally dozens of chemical confounds have been suggested at one time or 
.another as the substance of the Venus clouds. It was once assumed that the 
clouds were composed of water droplets, just as terrestrial clouds are, un- 
til spectroscopic studies revealed that the atmosphere of Venus contains 
very little water vapor. Spectroscopic data on the relative abundances of 
CO 2 , CO, HCJ,, HF and H 2 O in the atmosphere, as well as upper limits esti- 
mated for other proposed cloud constituents, have narrowed down the field 
of possible Venus cloud substances to several strongly hygroscopic materials 
such as carbon suboxide (CaOa) , an aqueous solution of hydrochloric acid 
(HCX..nH20) or an aqueous solution of sulfuric acid {H 2 S 04 .nH 20 ) • Accord- 
ing to Hansen and Hovenier (1974) , the best fit among these three proposed 
constituents to their cloud-droplet refractive indices deduced from polar- 
ization measurements is that given by a 75 to 83 percent solution of sul- 
furic acid, although within their given error limits a strong (about 30 
percent) solution of HC£ cannot be completely ruled out. As A- Young 
(1975a) has pointed out, experimental information about strong sulfuric acid 
or hydrochloric acid aqueous solutions, possibly contaminated by HF, at 
conditions similar to those in the Venus clouds, is sketchy or nonexistent. 
Therefore the cloud-composition problem may not be completely settled for 
a number of years. 

3.2 Atmospheric Motions 

Measurements of Venus winds have been made directly only by the 
Soviet probes Venera 4 and Veneras 7 through 10 during their descents through 
the atmosphere. UnfortTonately these measurements were made only up to heights 
of 50 to 60 km above the surface. Above this region there exist two main 
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types of wind measurements: (1) tracking of light and dark features appar- 
ent in ultraviolet photographs taken at low resolution from Earth or at 
high resolution from Mariner 10, (2) estimation of wind velocities from 
Doppler shifts of spectral lines. Figure 12 shows the results of all types 
of observations of Venus atmospheric motions to date. 

Veneras 4 and 7 landed on the dark side of Venus near the morning 
terminator, while Venera 8 landed on the sunlit side of the same termina- 
tor. Velocity profiles obtained by these three probes are quite similar, 
with near-zero zonal winds near the surface, gradually increasing with 
height up to 35-40 km, then strongly increasing to maxima of -110 to -150 

— T. *■ " 

m sec at a height of 50 km, where the measurements end. [A negative 
zonal wind velocity for Venus indicates retrograde motion, i.e., in the 
same direction as the planet's rotation.] 

In Figure 12 the regions of strong turbulence represented by shaded 
areas were found in radio occulation experiments performed with Mariners 
5 and 10 (Woo, 1975) . In both these experiments the path of the radio beam 
passed through the local late-night and early morning regions of the Venus 
atmosphere. This orientation corresponded quite well with the landing 
sites of Veneras 4 through 8. The lower turbulent region coincides with 
the strong wind shear measured by the Veneras; this suggests that the strong 
wind shear near 50 km is the cause of the turbulence in this lower layer. 

It is interesting to note that the wind profiles measured by 
Veneras 9 and 10, whose landing sites were nearer to the subsolar region, 
do not show evidence of the high velocities and strong wind she^s measured 
by Veneras 4 through 8 near the morning terminator. 

. The earliest evidence for Venus atmospheric motions came from earth- 
based telescopic observations in ultraviolet light (Boyer and Camichel, 
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1961; Boyer and Guerin, 1969; Dollfus, 1975). The prominent large, dark 
Y-shaped feature at the planet's equator appeared to have a remarkably con- 
sistent rotation period of about 4 days in a retrograde direction. Before 
radar observations established the true rotation period of Venus as 243 
days, this 4-day period had been assumed to represent the actual rotation 
of the planet. Since this was not the case, the only other possible con- 
clusion was that the whole Venus atmosphere at the cloud-top level near 
the equator was moving with an average velocity of -110 m sec The 
range of wind velocities which have been estimated from earth-based W 
photographs of Venus is indicated by the bracketed area in Figure 12. The 
double arrow labeled "4-day rotation" indicates the average velocity of 
the Y-shaped dark feature - 

Taking into account these apparent high wind velocities in the up- 
per part of the Venus atmosphere, Ainsworth and Herman (1977) have extra- 
polated the Venera 7 3^ situ measurements into two possible upper-atmosphere 
wind profiles, indicated by the dot-dashed cujrves in Figure 12. Profile 
(1) would be consistent with the suggestion by A. Young (1975b) that the 
estimated -100 m sec ^ winds in the Venus upper atmosphere are the re- 
sult of observational bias and that the actual winds in this region are of 
the order of no more than 5 m sec Profile (2) assmes that the strato- 
spheric wind velocities estimated from UV photographs and spectroscopic 
studies are not spurious results; in this profile the maximum velocity of 
-130 m sec ^ measiired by Venera 7 at 50 km is expected to persist well into 
the stratospheric region. 

High-resolution ultraviolet images from Mariner 10 support the gen- 
results of earth-hased UV studies. Belton ^ al . (1976a,b) were able 
to identify many details in these images in addition to the large dark Y. 
These smaller features included bow-like waves, circumequatorial belts, 
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a polar ring, and cellular features with lifetimes smaller than 4 days. In 
his vejy thorough analysis of the Mariner 10 images Limaye(1977) tracked 
many of these UV features having spatial extents of the order of 100 km over 
a period of two and one-half days. He estimated that these features are 
located at altitudes of 65 to 78 km above the planet's surface, and that this 
is a region of great atmospheric stability, i.e., the stratosphere of Venus. 
From a large sample of wind vectors distributed over most of the planet's 
area, Limaye found a mean zonal velocity component of -96.0 ± 16.0 m sec ^ 
and a mean poleward meridional component of 3.1 ± 10.2 m sec He also 
found that the zonal winds generally slowly increase in magnitude away from 
the equator to about 30 degrees latitude, then decrease rapidly toward the 
poles. The meridional profile of the zonal component of velocity can be de- 
scribed by the ecyiation 

u{({») = {sec (j> (1 - sin^fji) + cos <J) sin*^(|i} (3.1) 

where u{(j)) is the latitude-dependent zonal component of the wind, is the 

equatorial zonal wind, (p is the latitude, and a is an empirical constant. 

When a = 0, u(i)) = U cos (j), indicating solid-body rotation of the atmosphere 

o 

which requires that the angular velocity be constant with latitude. At the 
other extreme, in the limit as C£ approaches infinity, u((j)) == ssc cj> indicat 
ing that angular momentum is constant with latitude. Limaye has empirically 
determined that the value of a most consistent with his observations is 1.6, 
which means that the upper atmosphere of Venus is nearly in solid-body rota- 
tion. For this value of a, maximum zonal velocities of -101.1 m sec ^ occur 
at about 28° latitude, assuming a meridional distribution symmetric in lati- 
tude with an equatorial zonal velocity of .-95.1 m sec The asstamption of 
symmetry neglects the fact the zonal velocities measured for the southern 
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hemisphere were slightly higher than those for the northern hemisphere. How- 
ever/ the significance of these hemispheric dxfferences was lost in the un- 
certainty of the data. 

The result that the Venus stratosphere appears not to be exactly in 
solid-body rotation indicates that there must be a small non-zero meridional 
component of motion. The uncertainties involved in the tracking of the DV 
features, particularly regarding the true location of the poles, produce very 
large relative errors in the estimates of the meridional component of motion- 
The resulting standard deviation is several times larger than the average esti- 
mated meridional velocity, which' is small and poleward-directed - 

Limaye also analyzed the longitudinal variation of the zonal veloc- 
ities, particularly in the neighborhood of the sub-solar region. He found 
a general trend that zonal velocities between the sunrise terminator and the 
sub- solar point (local "morning" on Venus) tended to have smaller than average 
magnitudes, while velocities measured for local "afternoon" regions were 
larger than the average. The velocity increase from morning to afternoon 
was of the order of 10 m sec ^ . Given the retrograde motion of the Venus 
atmosphere, it is obvious that features eastward (toward sunset) of the sxob- 
solar point have been in the sunlit hemisphere of Venus longer than those 
nearer to the morning terminator; this extended opportunity for solar heating 
correlates with the increase of zonal velocities during the day. 

An increase in average zonal velocity from morning to afternoon on 
Venus was also noted by Traub and Carleton (1975). These investigators 
made spectroscopic measurements of wind velocities using a CO 2 absorption 
line at 8708 S and a Fraunhofer line at 6716 S. The zonal wind velocities 
were determined by comparison of Doppler shifts of the same spectral line 
measured while viewing two different areas of the illuminated portion of 
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the planetary disk. The areas viewed were usually located on the equator 
at either the edge of the bright limb or near the terminator. Zonal veloc- 
ities estiniated in this manner for the CO 2 line varied from -2 to -125 m 

-1 -1 

sec , wxth an overall mean value of -83 m sec , However, averaging local 

morning and afternoon measurements separately, they found generally smaller 

. . . -1 -1 
morning velocities (-73 m sec ) than afternoon velocities {-111 m sec ) . 

These results are indicated by the points labeled Tcm and Tea in Figure 
12. For the Fraunhofer line, the results represent an average of only 
four measurements made viewing the morning region of Venus. The average 
velocity obtained (-78 m sec refers to the motion of the opaque cloud 
deck from which the solar photons are reflected. The level of CO 2 line 
formation (for the narrow core of the line used for determining its precise 
position) appears to be somewhat above the main cloud, and the morning 
velocity determined from the CO 2 line is about 10 m sec ^ slower than that 
estimated for the main cloud. This appears to be a marginally significant 
result suggesting that in the morning hemisphere zonal velocities are in- 
creasing with depth into the Venus atmosphere. 

Preliminary results from independent spectroscopic measurements of 
Venus zonal wind velocities have been very recently reported by Schorn and 
Young (1977) . Their measurements of Venus spectra in the 7600 £-8800 £ 
region have yielded a mean rotational velocity of 40 ± 47 m sec ^ retro- 
grade, which is considerably smaller than the mean values reported by Traub 
and Carleton, but not outside the range of their individual measurements. 

Theoretical models of the dynamics of the Venus atmosphere have had 
some difficulty explaining the observed temperature structure and motions 
of the entire atmosphere. Stone (1975) has examined the time scales char- 
acteristic of important atmospheric processes on Venus. His results. 
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together with observed wind velocities and temperatures, indicate that the 
Venus atmosphere contains two distinct regimes. The deep atmosphere, below 
40 km, is characterized by a near-adiabatic taaperature lapse rate, weak 
zonal , meridionals and vertical velocities , and negligible diurnal effects . 
The more successful dynamical models assume that the greenhouse effect in 
this region is sufficient to cause the large observed surface temperatures 
and that the atmospheric motions up to about 40 km are due to a Hadley- 
type circulation (Sasamori, 1971; Kalnay de Rivas, 1975; Gierasch, 1975). 
These models can adequately explain the adiabatic, nonturbulent stsnicture, 
the lack of strong horizontal gradients and the observed small horizontal 
velocities in the troposphere of Venus. 

In the upper atmosphere, above about 70 km a different dynamic 
regime predominates; diurnal effects appear to be important, the temperature 
lapse rate is sub- adiabatic, and the zonal motions are quite strong. Ac- 
cording to scaling of the equations of motion carried out by Leovy (1973) 
and others, the large-scale motions in the upper atmosphere appear to be 
in cyclostrophxc balance. This can also be seen in mosaics of pictures 
from idariner 10; the main circulation pattern in each hemisphere is a cir- 
cumpolar vortex (Limaye, 1977). However, no mechanism yet proposed satis- 
factorily explains how this circumpolar vortex is formed in the Venus stra- 
tosphere. One likely possibility is the ’’moving flame" mechanism first 
studied in the laboratory by Fultz ^ (1959). The moving flame, which 

jLn the case of the Venus stratosphere would be strongly localized heating 
near the bottom of the layer due to CO 2 absorption in the sub-solar region, 
induces a thermal wave in the fluid. This wave lags behind the moving heat 
source by, an amount which increases with height, because of the finite 
time required for conduction to transmit the heat into the fluid. There 
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.soon arise convection cells which also have a tilt increasing with height 
from the bottom of the layer. Because of this tilt, horizontal momentum 
in a direction opposite to that of the motion of the heat source is trans- 
_ported upward. Thus the net motion at the top of the fluid layer has a di- 
rection opposite to that of the motion of the localized heat source, which 
-in the case of Venus is the sun. 

_ . . Numerical computations by Young and Schubert (1973) seem to indi- 

cate that the "movxng flame" mechanism alone can account for both the 
magnitude and direction of the observed 4-day rotation of the Venus upper 
atmosphere. Essential. to this result is the assumption that the atmosphere 
is stably stratified above altitudes of about 60 km, which is apparently 
the case in the real Venus stratosphere. However, as Stone (1975) points 
out. Young and Schubert were not able to obtain numerical results for the 
large values of the parameter A (the square of the inverse Froude number) 
•appropriate for the real Venus atmosphere, since their iterative technique 
failed to converge for A values an order of magnitude smaller than that for 
Venus. Extrapolation of the results they did obtain for smaller values of 
A does produce zonal velocities within a factor of 2 of the observed wind 
speeds. Nevertheless, this extrapolation is by no means conclusive evi- 
dence that the "moving flame" is responsible for the Venus- stratospheric 
circulation. 

iinother mechanism which may be partially responsible for the Venus 
upper atmospheric circulation has been discussed by Fels and landzen (1974) . 
Internal gravity waves generated by diurnal heating will tend to transport 
prograde momentum upward out of the stratospheric layers where the absorp- 
tion of solar radiation occurs, leaving these absorbing layers with net 
retrograde momenttim. Prograde velocities supported by this mechanism cannot 
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exceed the apparent speed of the sun, about 4 m sec . Retrograde veloc- 
ities of the magnitudes observed, about 100 m sec can be supported only 

-X 

if a mean retrograde velocity of about 25 m sec is produced independently 
by some other mechanism, filso, these strong retrograde velocities would 
exist only in a layer at most 5 km thick. The two-dimensional model of 
Eamanathan and- Cess (1975), which neglected rotational effects, also 
showed that propagating internal gravity wave generated by diurnal heating 
could induce mean zonal velocities in the Venus stratosphere. Their anal- 
ysis differs from that of Pels and Lindzen in that the diurnal heating oc- 
curs only in the upper part of this atmospheric region above altitudes of 
about 85 km. This diurnal heating and the associated gravxty waves would 
induce mean zonal velocities in both the lower and upper parts of the stra- 
tosphere. Thxs model results in retrograde mean zonal velocities of the 

order of those observed in the lower stratosphere, with these mean retro- 

‘ - -1 
grade velocities increasing with height to about 300 m sec near 100 km. 

However, the model results also include a ditumally fluctuating velocity 
component above' 85 km which may have either prograde or retrograde direc- 
tion and a magnitude exceeding that of the mean velocity. Thus, depending 
on the local time of day on Venus, the net velocities in the upper strato- • 
sphere can be prograde in direction. 

More recently Young and Pollack (1977) have developed a three- 
dimensional model for the Venus atmospheric circulation from the planet's 
surface up to about 65 km. The principal physical mechanisms active in 
their model appear to be the mean merxdio’nal circulation, which initially 
transports net retrograde momentum upward from the slowly rotating plane- 
tary surface, and a nonlinear instability, which may be of finite amplitude, 
involving both the mean meridional circulation and planetary-scale eddies 
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arising from dinrnal temperature variations in the upper atmosphere. The 
model was able to reproduce in direction and magnitude the observed mean 
zoneil wind velocities of about -100 m sec ^ in the stratosphere, as well 

as other obseirved properties of the circulation such as small horizontal 

. ' ~1 
temperature gradients and meridional velocities of the order of 10 m sec 

or less in low and mid-latitudes. Under certain conditions the model also 
produced planetary waves which exhibited vertical flow field contours in 
a horizontal Y-shaped configuration. However, although the agreement with 
the observed circulation appears to be good, the results of this model 
cannot be considered to be conclusive, particularly in regard to the vari- 
ability of the Venus zonal winds with time and location, as observed by 
Traub and Carleton (1975) and Limaye (1977) • The vertical profile of mean 
zonal velocity (averaged over all longitudes) deteimiined from the model 
appears to be very similar to the Venera 8 profile (Figure 12) . However, 
that profile refers only to a single region of the atmosphere near the 
morning terminator and is not necessarily representative of the veolcity 
profiles for other longitudes on the planet^ For example, the Venera 9 
and 10 velocity profiles, which refer to longitudes near the sub-solar 
point, appear to differ significantly from that measured by Venera 8, es- 
pecially in the upper part of the atmosphere where diurnal effects appear 
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to be more important. It is hoped that the planned modifications in the 
Young and Pollack model to include more realistic treatment of radiative 
heating and cooling and to extend the upper boundary to about 80 km will 
result in more realistic velocity profiles at all longitudes and will pro- 
vide more insights into the general circulation in the atmosphere well above 


the clouds. 



CHAPTER IV 


ESTIMATION OF VENUS ZONAL VELOCITY PROFILES 
USING PRESSUEE-INDUCED LINESHIFTS 

4»I Description of the Spectra 

As mentioned in the introduction to this paper, the Venus spectra 
analyzed in this work were measured by Dr. P. Connes using an improved 
version of the Connes’ interferometer at the coude’ focus of the Hale 5 m 
(200 in) telescope at Mt. Palomar in October 1973 and June 1974. The best 
of these measured spectra were incorporated into two averaged spectra 
identified as Venus A and Venus 5. 

Venus A is the average of four recordings made October 14, 15, 22 
and 25 in 1973. Each recording was made over a two- to three-hour period 
beginning at about 2:00 P.M. local time. The planet appeared as an even- 
ing star with an average phase angle (defined as the Venus-centered angle 
between the sun and the earth) of 11 . During the period of the observa- 
tions the disk of Venus had an average apparent semidiameter of 9.8 arc 
sec. 

The spectrum identified as Venus 5 is the average of five record- 
ings made on consecutive days from June 12 through 16 in 1974. Each of 
these spectra was also recorded over a two- to three-hour period, beginning 
at approximately 4:00 A.M. local time, since Venus was appearing as a morn- 
ing star. Its average phase angle was 56.6°, and its average semidiameter 
was 6.9 arc sec. Specific parameters for all the individual recordings are 
given in Table 3. In each case the field of view was a circular area 12 
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arc sec in diameter centered on the midpoint of the illuminated limb. A 
more detailed discussion of the viewing geometry is given in Section 4.2. 

For calibration purposes, during the measurement of the spectra 
making up Venus 5, an absorption cell containing a known amo\int of CO at 
low pressure was placed in the path of the interferometer. This caused 
sharp, unshifted CO absorption lines to appear in the Venus 5 spectirum 
along with the Doppler- shifted Venus CO lines. Lines in the 2-0 band of 
C^^O^® due to the abso 2 :ption cell were used to calibrate the Venus 5 spec- 
trum, using as st^dards vacutim line positions for tliis same band measured 
with a similar interferometer in the laboratory by Guelachvili (1973) . The 
positions of 23 telluric absorption lines due to (0-0) electronic transi- 
tions of Qz near 7850 cm ^ were determined from the calibrated Venus 5 
spectrum; the positions of these Oj lines were then used as standards to 
calibrate the Venus A spectrum. After calibration, line positions for un- 

blended telluric O 2 lines in this spectral region (7820-7900 cm ) measured 

—3 —1 

in the two averaged spectra agreed within ± 0.4 x 10 cm 

A preliminary survey of the spectra was carried out to determine 
whether pres sure- induced lineshifts in the HF 1-0 and 2-0 bands and in the 
HCA 2-0 band due to the CO 2 in the Venus atmosphere could actually be de- 
tetected. Absolute wavenumbers for these three absorption bands measured 
by Guelachvili (1976) with a similar Fourier-transform interferometer were 
used as standards for comparison to the Connes planetary spectra- The stan- 
dard line positions were corrected for Doppler shifts due to the relative 
motions of Venus and Earth and also due to the rotation of the earth. These 
corrected standard line positions were then compared with the HF and HCil 
line-center positions measured in the Venus A and Venus 5 spectra- The 
small residual lineshifts were then compared with laboratory measurements 
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of lineshifts due to CO 2 (Guelachvili and Smith, 1978) ; these results are 
shown in Figures 15 and 16. 

In the two figures measurements belonging to poor-quality Venus 
spectral lines are indicated by open symbols, and "good" lines are denoted 
by solid symbols. All the HF and HOI lines identified- in the two averaged 
spectra are listed in Table 4 along with their "quality factors"; a quality 
factor of 5 indicates that the spectral line is unblended and almost to- 
tally free of noise, while a line having a zero quality factor is identi- 
fiable in the spectrum but totally unusable for estimating pressure-induced 
lineshifts. The relatively large number of poor-quality lines is due to 
strong telluric absorption in some spectral regions such as the P-branch 
of the HF 1-0 band, and also to blending of several of the hydrogen halide 
absorption lines with some of the many Venus CO and CO 2 lines in the spec- 
tra. The net resiilt is that both of the HF bands contain so few good lines 
that they cannot be used for the study of atmospheric-xnduced lineshifts; 
the P-branch of the HC£ 2-0 band is also of poor general quality. How- 
ever, in the R-branch of this band, nearly all the lines belonging to both 
isotopes are quite usable, having quality factors of 3 or greater. A por- 
tion of each spectrum containing several of the best-quality lines in the 
HCX- band is shown in Figure 14. 

J 

Even the uncertain data points can provide some information, and 
comparison of all the measured residual lineshifts with the pressure-induced 
shifts by CO 2 measured in the laboratory indicates that pressure- induced 
lineshifts of HF and HcjZ, perturbed by the CO 2 in -the Venus atmosphere are 
large enough to be detected in these high-resolution planetary spectra. 

This is best seen in the R-branch of the HCfi. 2-0 band, where the Venus 
spectra show residual lineshifts increasing with J in the same manner as 
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the laboratory measurements of the pressure~induced lineshifts. Moreover/ 
the entire Venus lineshift cu 2 nres are apparently displaced by a constant 
amount from the laboratory cur;ye; this indicates a small residual Doppler 
shift which is apparently due to motions ^n Jbhe Venus atmosphere, since 
the planet's true rotational velocity is nearly zero. 

- With the aid of a computer program called DEEPTH developed by re^* 
searchers at the Centre National de la Recherche Scientif ique , Orsay, it 
is also possible to study in detail the symmetry of each absorption line. 

The program calculates the central positions of chords drawn across the 
line at five levels equally spaced from the base to toe peak of the line 
(Figure 17) . For the HCA lines studied in detail, it was found that the 
centers of. these chords were increasingly displaced from toe peak to the 
base of each line. Telluric lines in thxs same spectral region do not show 
the same asymmetric shape as the HCil lines. Considering these asymmetries 
together for the whole R-branch of the 2-0 band, the displacements of toe 
chord centers appear to represent not only shifts due to increased CO 2 pres- 
sure at lower levels in the Venus atmosphere , but also a residual 'Doppler 
shift due to the variation of Venus wind velocities with height. The par- 
ticular nature of the shift pattern for HCJl makes it possible to calculate 
profiles of Venus zonal wind velocity as a function of atmospheric pressure, 
using the line asymmetries measured in the R-branch of the HC^. 2-0 band. 
Details of the method of calculation and the resulting velocity profiles 
will be discussed in later sections of this chapter. 

4.2 Geometry of the Problem 

In order to estimate wind velocities on Venus from small residual 
lineshifts, it is first necessary to calculate as accurately as possible 
the contributions to the lineshifts which arise from the motions of Earth 
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and Venus in their respective revolutions about the Sun and from the rota- 

- tion of the Earth. Also, since the spectra contain contributions from a 

- rather large area of the illuminated portion- of the Venus disk, we must 

- take into accovmt the variation of intensity of illumination of the planet 

- and the latitudinal variation of zonal velocities on Venus in order to trans~ 
late observed residual Doppler shifts of the spectral lines into zonal 

wind velocities in the planetary atmosphere. 

4.2.1 Effect of translational and 
rotational Doppler shifts 

The largest contribution to the shifts of the Venus spectral IdLnes 
comes from the relative orbital motions of Earth and Venus. These trans- 
-lational Doppler shifts have been calcinated and tabulated for convenient 
reference by Niehaus and Petrie (1961) ; Figure 18 is a graphical represen- 
tation of their calculations for Venus for the period 1971 to 1977, It 
can be seen from the figure and Table 5 that the translational Doppler 
shift remained fairly constant throughout the 11-day period of observations 
for Venus A, while for the Venus 5 observations the shxft factor decreased 
more rapidly in magnitude over a period of 5 days. The net change in the 
translational Doppler shift factor from the first to the last day of ob- 
servations is about +1.2 percent for Venus A and about -2-1 percent for 
Venus 5. The larger variation in the translational Doppler shift causes 
the same lines in the Venus 5 averaged spectrum to be somewhat broader 
than those in the Venus A spectrum- 

Since the Niehaus and Petrie (1961) calculations of the transla- 
tional Doppler shift were intended to be a rough estimate for quick refer- 
ence during observations and are accurate to only two signifxcant figures, 
it is necessary for detailed study of the spectra to calculate the 
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translational shift more acrnrataly. These calculations were made using 

listings from the American Ephemeris and Nautical almanac (1973, 1974) 

giving the true distance between Earth and Venus at 0 h Universal Time 

for each day of the observations. The Doppler shift factors, listed in 

Table 5, were calculated for the midpoint of the two- to three-hour time 

period during which each interferogram was taken. 

A smaller component of the Doppler shift, which arises from the 

Earth’s rotation, depends on the specific location of the observatoiry and 

on the average position of Venus in the sky at the time of each observation. 

As listed in. Table 1, the tangential speed of a point on the earth's sur- 

-1 

face at latitude ^ in m sec is 463.8 cos <j». Since Mt. Palomar is located 
o 

at 33.18 N (geocentric latitude), its tangential velocity is about 387.9 

m sec Thus the maximum possible rotational Doppler shift factor, if 

Venus were located at the horizon due east of the observatory, would be 
-5 

+0.129 X 10 , which is an order of magnitude smaller than the transla- 

tional Doppler shift. 

The component of the earth's tangential velocity in the direction 
of Venus can be expressed as 

v^ = -U^ cos a sin z (4.1) 

where U^ is the eastward tangential velocity of the observatory, and a and 
z are respectively the average azimuth (from due north) and zenith angles 
of the planet Venus as seen from the observatory during the measurement of 
the interf erograms . Using spherical trigonometry, a and z can be calculated 
as functions of the latitude and longitude of the observatory, the time of 
day the observations were made, and the apparent position of Venus in the 
sky at 0 h Universal Time as given in the Ephemeris for each observation 
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date. The results of these calculations for the individual Venus spectra 

4 

are listed in Table 5, along with the net Doppler shift factors combining 
both translation and rotational effects. 

4.2.2 Scaling of residual Doppler shifts 
to obtain Venus equatorial zonal winds 

In order to calculate Venus equatorial zonal wind velocities from 
the small residual Doppler shifts measured in the spectra, it is necessary 
to take into consideration the spherical shape of the planet, the varia- 
tion of the zonal velocities with latitude, and the variation of bright- 
ness over the sampled portion of the planetary disk. The average areas 
viewed for each set of interferograms are shown in Figure 19. From the 
residual Doppler shifts determined from the spectra we can readily calculate 
a velocity v^, which is the eeirthward component of the Venus atmosphere's 
apparent rotational velocity averaged over the obseirved circular area of 
the illuminated part of the planet's projected disk. Since at no time 
during the observations does Venus appear in full phase (when the sub- 
earth and sub-solar points on the planet coincide) , and the viewing aper- 
ture is not centered on the sub-earth point, the usual air-mass parameter- 
izations which describe the relative contribution of light reflected from 
each point on the planet to tdie spectrum (Moroz, 1967) , are not very help- 
ful in this case. 

The average earthward velocity component v^ can be expressed as a 
function of the planetary geometry, zonal velocity distribution and bright- 
ness distribution by the equation 
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(4.2) 


where u(({!) is the zonal velocity distribution as a function of latitude, 
I„(A,(j)) is the brightness distribution relative to that at the sub-solar 
point, and the sin A term arises from the projection of the spherical 
planet onto the plane of the apparent disk. The variables A and (f) are re- 
spectively longitude and latitude on Venus, with the origin defined as the 
sub-earth point, i.e. , the center of the apparent disk- 

The simplest velocity distribution which can be used on the above 
equation is that corresponding to uniform solid-body rotation of the atmo- 
sphere, i.e., u(<ji) = U cos <|), where U is the equatorial zonal velocity. 

o o 

However, since the integration of (4.2) is done numerically, it is not very 
difficult to use a more realistic expression for the velocity distribu- 
tion such as that determined empirically by Limaye (1977) , which is given 
in Equation (3.1) . The zonal velocity can be expressed as the product of 
the constant equatorial velocity with a zonal velocity distribution 

fvinction u(<^). The constant U can be taken out of the integral in the 

o 

numerator, allowing the calculation of a convenient velocity scaling fac- 
tor = U^. Velocities obtained from residual Doppler shifts in a spec- 
trum can then be quickly divided by the appropriate scaling factor to 
obtain the average equatorial zonal wind velocities. 

Several distributions for the solar illumination of the disk of 
Venus were considered, but for the sake of simplicity of calculation, only 
Lambert's Law for diffuse reflection from a solid surface and the reflection 
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law for isotropic scattering seemed appropriate for use in the integration 
of (4.2). Exact calculations of the brightness distribution for Venus as- 

y * r , 1 . ' -z. 

suming a semi-infinite cloud with Mie scattering have been carried out 
using the method of- spherical harmonics by Devaux et al. (1975) for wave- 
lengths ranging ifemT the visible orange (A = 0.6 pm) to the ultraviolet 
(A = 0.4 pm). However, Hapke (1976) has pointed out that calculations as- 
suming ‘isotropic 'scattering by 'the— Venus- clouds-are an~reasonable"agree- 
ment with observations for this' same wavelength region, particularly to- 
Wcird the orange wavelengths, -where -the Tsi'ngle- scattering 'albedo of the 
clouds approaches 1.00. In the infrared region of interest in the present 
work (1.5 to 2.0 pm), -the ’ single-scattering albedo is also close to 1.00, 

^d the isotropic- scattering 'regime is applicable. A mathematically simpler 
distribution, which may be a good first approximation to -isotropic scatter- 
ing, is that given 'by Lambert ' s 'Law for 'diffuse reflection, which assumes 
that the reflected light 'as isotropic and-is independent-of-the state of 
polarization and the angle -of incidence of the incxdent light. The bright- 
ness-distributions due to both 'tjpes'of 'diffuse reflection 'laws .will be com- 
pared here for the Earth- Venus-Sun- configurations at' the time of the measure- 
ments for Venus A and Venus 5. " 

■ * ' Following the notation used by Chandrasekhar (1960), 'the intensity 

of illumination I(p,p^) of any point M on the planet’s surface as seen from 
space is the function of the cosines of the zenith angles of the observer 
(p) and the sun (p^) as seen from point M. These zenith angle cosines in 
turn can be expressed as functions of the Venus phase angle (V^) and of the 
latitude and longitude of point M, that is, 

p = cos A cos (4.3a) 


pQ= cos (A-V^) cos (j> 


(4.3b) 
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vrfiere as before, longitude and latitude (X,cj)) are defined with the sub- 
earth point as the origin, and is taken to be positive when Venus is 
an evening star and negative when it is a morning star. Lambert's Law 
for diffuse reflection can then be stated as 


) = w u p 
o o o 


(4.4) 


where is the single- scattering albedo and F is the incident flux at 

the top of the atmosphere at the sub-solar point, where p =1. Thus 

o 

the brightness of any point on Venus relative to that of the sub-solar 
point is simply the cosine of the solar zenith angle at that point, i.e,. 
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(4.5) 


The diffuse reflection law for isotropic scattering, however, is 
more complex, and the brightness at any point on the planet depends on the 
positions of both the earth and the sun in relation to that point. The 
law is given by 


0) F vt 

I(U,y^) = — ( 

4 y+y^ 


) H(y) H(y ) 
o 


(4. '6) 


where H(y) and H(y^) are Chandrasekhar H-functions which also depend on the 

value of the single-scattering albedo (i) . The values of the H-functions 

o 

for given values of y and can be found by lengthy calculation or in al- 
ready published tables (Chandrasekhar, 1960). 

Using the tables of H-functions, the longitudinal distribution of 
relative brightness of the disk of Venus along the equator and at two other 
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latitudes,, assuming isotropic scattering, was calculated for both the Venus 
A: and Venus 5 viewing configurations. Figure 20 shows the calculated 
brightness- distribution for isotropic. scattering (dashed lines) , as well 
as the. corresponding calculated distribution assuming lambert's Law re- 
flection (solid lines) . It is apparent in the figure that the results for 
the two laws' of diffuse reflection are nearly identical, except near the 
extreme edge of the planetary disk in lower latitudes. Therefore it seems 
reasonable to use the mathematically simpler Lambert's Law to describe 
the brightness distribution of the planetary disk for the particular view- 
ing configurations of the spectra studied in this work. Since most of the 
differences between the two distributions appear near the sub-solar region, 
it is expected that the more realistic isotropic- scattering distribution 
should be used when Venus appears at or near full phase. 

It should also- be noted that in Figure 20 the coordinates given for 
location on the planet are not latitude and longitude, but Cartesian co- 
ordinates in the plane of the apparent disk' of Venus with the origin at the 
midpoint of the bright limb; horizontal and vertical coordinates (x,y) are 
scaled so that the planetary radius S = 1. The sub-earth point in this 

system then has the coordinates (-x ,0) . This transformation to Cartesian 

o 

coordinates facilitates the integration of the velocity scaling function 
over the illuminated portion of the circular field of view. Careful com- 
parison of the Venus phase angle cind semidiameter for each observation date 
reveals that for a number of the individual spectra the field of view, which 
is circular with a more or less constant diameter of 12 arc sec, overlaps 
both the terminator and the edge of the planetary disk (Figure 18b) . There- 
fore the area over which integration of Equation (4'. 2) is performed may not 
be a circle, but a circle with arcs of larger radius cutting off its edges. 



39 


Figure 21 shows the distribution of the relative contribution to the spec- 
troscopically observed earthward velocity component for both Venus A (a) 
and Venus 5 (b) , isoplethed for the observed illtuninated area of the disk. 
'The quantity isoplethed in the figure represents- the integrand in the 
numerator of the following equation, which was evaluated numerically 


V 

C 

F =5 — = 

V U, 


//u(y) (x+Xq) I^(x,y) dx dy 


(4.7) 


in order to find the velocity scaling factor F^ for each of the individual 

spectra and for the averaged spectra. In the equation above u(y) , (x+x ) 

o 

and I (x,y) are the Cartesian-coordinate equivalents of u {({)), sin X and 
Ij^(X,<{)) used in Equation (4.2). The results of the numerical integrations 
show that the velocity scaling factor F^ is 1.147 for Venus A and -0.873 
for Venus 5. The scaling factors for the individual spectra differ from 
these mean values by not more than ± 0.029. 


4.3 The Direct Problem: Calculation of 

Expected Line Asymmetries for a Model 
Atmosphere 


Direct calculation of expected line shapes in the R-branch of the 


2-0 band for given profiles of temperature, pressure and equatorial 
zonal wind velocity (U^) in the Venus upper atmosphere aids in the under- 
standing of the relationship of the pressure- induced and wind-induced line- 
shifts at specific atmospheric levels to the asymmetric profile of the en- 
tire absorption line. Such a calculation is relatively straightforward, 
using the standard formulas for the absorption coefficient and transmittance 
for a line having the Voigt profile in a plane-parallel, non-homogeneous 
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atmosphere (Goody, 1964)-. In this model it is assumed that no significant 
scattering of the infrared radiation occurs in the atmosphere above the 
clouds, and that the cloud tops act as an unbroken surface reflecting in- 
frared radiation according to Lambert's Law. 

The variation of intensity for an incident ray of intensity 
passing through a region in thermodynamic equilibrium containing an amount 
of absorber dm is given by 

where k^ is the coefficient of absorption and is the black-body inten- 
sity at the temperature of the medium. In the case of sunlight reflected 
from the Venus cloud tops, I^ is at least 10 orders of magnitude greater 
than Thus we can neglect B^ in the transfer Equation (4.8) and inte- 

grate to obtain 


I = 



e 



(4.9) 


where I^ is the original intensity of the ray before it has passed through 
o 

the medium. We can also define the transmittance as = 1^/X^ = 

o 

exp C - / k dm) . 

V • 

The absorption coefficient k^ depends not only on frequency, but 
also on the pressure, temperature and molecular composition of the absorb- 
ing medium. At the conditions present in the Venus atmosphere 50 to 80 
km above the surface, the HC^ absorption lines are expected to have a 
Voigt profile. Calculation of expected magnitudes of both Doppler-broadened 
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and collision-broadened half-widths , as well as line intensities and mean 
amounts of HC5, absorber for this region of the Venus atmosphere shows that 
the HCA absorption lines should fall in the linear and transition regions 
of the curve of growth (see Goody, 1964, p. 237^ . The Voigt line shape can 
be expressed as 



dy 


(4.10) 


where S is the line intensity, is the Doppler half-width at 1/e height 
of the line, and a = a / a , where a is the Lorentz, or collision broadened 

Xj D Xj 

half-width at half-height of the line. Under the integral y is a dummy 

variable of integration, while t = (V-V )/ a„, where V is the central fre- 

o D o 

quency of the line. 

The value of a depends on both pressure and temperature as follows 

Xf I 


ct. 




(4.11) 


where a is the Lorentz half-width measured at standard pressure p (usually 
s s 

1 atm) and standard temperature (usually about 300 K) . The exponent n 
is generally taken to be 1/2, which is the value obtained assuming that molec- 
ular collisions occur with the frequency given by classical kinetic theory. 
Unfortunately, complete measurements of C02“broadened Lorentz half-widths 
in the HCS< 2-0 band are for the RO and RIO lines by Rank et al. (1960a) . 
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However, reasonable estimates of CO 2 •"broadened halfwidths for ■the other 
lines in ■this band can be made using existing measurements of self- 
broadened halfwidths in this same band by To^th et al . (1970) . By compar- 
ing measurements in the fundamental of HC5, of self-broadened (Benedict et 
al. , 1956) and CO 2 -broadened (Varanasi et al. , 1972) halfwid-ths, one can 

calculate broadening coefficients y(C 02 ) = a (C 02 )/ct (HCA) for each of the 

s s 

lines in ■this band. Careful examination of the results of Rank ^ al . 

(1963) for broadening by Kr and Xe in both the 1-0 and 2-0 bands of HCJl 
shows that the ratio of broadening coefficients for lines having identical 
J values in the two bands, Y(2-0)/ y(l-O) remains more or less constant 
for different perturbing molecules. Thus the CO 2 -broadened halfwidths in 
the HCS, 2-0 band can be calculated by multiplying the measured self-broadened 
halfwidth for the 2-0 band by the broadening coefficient for the 1-0 band 
determined from the Varanasi ^ al . (1972) measurements, and finally by 
the previously determined ratio of the .broadening coefficient y(2-0)/ y(l-O) . 
The results of such calcula-hions were used along wi-fch the two direct measure- 
ments by Rank et al. (1960a) in Equation (4.11) in the model. 

The Doppler halfwid-th is independent of pressure, depending only 
on -the ambient temperature and -the central frequency of the absorption line. 
It is calculated directly -using the formula 



where is the wavenumber of -the line center in cm ^ , c is the velocity of 

■s 

light, k is Boltzmann's constant' and m is the mass of a single HC5- molecule. 

The line intensity S also depends on temperature, which affects the 
population of the lower energy state involved in the absorption transition. 
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For a line in the HCi 2-0 band the intensity can be expressed as 


S(T) = S{T ) 
s 


W 


exp [ ( 


1 

T 


)] (4.13} 


T 


where again T is the standard temperature at which the reference lane in- 
s 

tensity S(T ) was measured. The function Q (T) is the rotational oartition 
s R 

function for the molecule; its exact values have been tabulated for a num- 
ber of atmospheric molecules , but it can also be approximated by 
Qj^(T) = kT/hcB for sufficiently large values of and T. As pointed out 

by Herzberg (1950) , this approximation appears reasonable for HCS- at tem- 
peratures near 300 K. Under the approximation the ratio Q_(T )/ Q„{T) ap- 

Iv S XV 

pearing in Equation (4.13) would be replaced by T /T. 

s 

In Equation (4.10) the central frequency V is usually taken to be 

o 

constant for a given line. However, in this study of line asymmetries, 
the central frequency varies with pressure, through both the pressure- 
induced lineshifts and the Doppler shift due to the wind velocity at each 
pressure level, i.e.. 


Vo(p) 


(1 + f. 


(P) ) 

o 


+ S 

s 



(4.14) 


where fjj(p) is the Doppler shift factor Av/V arising from the wind velocity 

specified for each pressure level, V is the standard vacuum, wavenumber for 

s 

o 

the line center, and 6^ is the laboratory-measured pressure-induced shift 
per atm of the perturbing gas. 
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Since a number of terms in k^(p,T) are directly dependent on pres- 
siire, it is more convenient to express the incremental path length dm in 
Equation (4.9) also in terms of atmospheric pressure. In a plane-parallel 
Venus atmosphere, the amount of HCA in each increment of the path is 

P 

^ = q ri dz (4.15) 

Ps 

where q is the volume mixing ratio of HCil to CO 2 , p is the actual density 
of the atmosphere at height z, and is the density of CO 2 , which makes up 
95-99 percent of the atmosphere, at STP. The effects of non-vertical paths 
of the beam through the atmosphere are taken into account by calculating an 
average airmass factor rj for the viewing config^zration. For a beam of sun- 
light reflected toward earth by a single point on the surface of Venus, 

n = sec 0 + sec 9, where 0 and 0 are respectively the zenith angles of the 

o o 

sun and the earth at that point. For the whole field of view for each spec- 
trum an average airmass factor, weighted according to the relative illunina— 
tion of each point of the planet, is calculated by the following formula 
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(x,y) dx dy 


(4.16) 


where u = cos 0 and u = cos 0 and the Cartesian coordinates (x,y) as well 
o o ^ 

as the integration limits, are the same as those defined in Section 4.2.2. 

In the airmass calculation for the Venus spectra a Lambert’s Law brightness 

distribution [I (x,y) = y 3 was used, and the integrations were carried out 
R o 
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numericcilly in the same manner as for the calculation of the velocity scal- 
ing factor in Section 4.2.2. The results of these calculations give ri = 
3.427 for Venus A, and n = 4.804 for Venus 5, _ ■ 

Assuming that the atmosphere is generally in hydrostatic equilibrium, 
i.e./ dp = -pgdz, then the above equation becomes 


~ <1 - 

dm = n dp 


(4.17) 


where g, of course, is the acceleration of gravity for Venus (885.2 cm 
sec ) . Thus the transmittance of a given layer of the Venus atmosphere in 
the neighborhood of each HC5. line becomes 


Pi 


= exp [ 


-/ 


<1 n 


P2 




(p,T) dp ] 


(4. '18) 


Finally, in order to approximate as nearly as possible the actual 
shapes of the lines as they would appear in planetary spectra, the profiles 
were convolved with the instrumental line shape (ILS) appropriate for the 
planetary spectra. Since all the Venus spectra were unapodized, the proper 
ILS is simply the classical Fourier Transform of the apodizing function 
P(5) =1, where P(6) is the weighting given to the interferogram as a func- 
tion of the path difference 6. The ILS thus is of the form (sxn x)/x, where 
X = 2710 M; 0 is the distance in cm ^ from the center of the instrumental line 
profile, and M is the maximum path difference in cm. For the Venus spectra 
used in this study, the ILS is quite narrow, having a half-width at half-height 
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of 0.006 cm and convolution of this ILS with the theoretical profile 
does not significantly alter the line shape. An excellent discussion, along 
with graphical representations, of several types of apodizing. functions 
used, in Fourier spectrometry can be found in the work of Norton and Beer 
(1976) . 

In the model calculations, S (T) , (T) , a^(p,T) and V^(p) for 

the RO through R8 lines in the 2-0 band were computed for 46 levels 

spaced 1 km apart between 50 and 95 km above the surface. The pressure and 
temperature profiles used (Figure 22) were from the Venus model atmosphere 
given by Mairov and Ryabov (1972) . Four wind velocity profiles were used 
and are shown in Fxgure 22. The profile designated ”A-H(1)" in the figure 
is a linear representation of the extrapolation by Ainsworth and Herman 
(1977) of the Venera 7 data to allow small zonal velocities in the upper 
atmosphere (see also Figure 12). Likewise, the line "A-H(2)” is their other 
extrapolation of the Venera 7 results, specifying a constant zonal veloc- 
ity of -130 m sec ^ throughout the upper atmosphere. The "Median” model 
velocity profile also shows a decrease to near-zero wind speeds at 90 km, 
but was drawn so as to include the measurements of Traub and Carleton (1975) 
and Limaye (1977), as well as those of Venera 7. A fourth model, called 
"A-H(2B)" has the same constant velocity profile as "A-H(2)", with the ex- 
ception that the zonal velocity is zero in the layer -between 70 and 72 km. 

For each level in the model the absorption coefficient k(V,p,T) was 
calculated at discrete frequencies spaced 0-002 cm I apart, extending over 
a region ± 0.100 cm ^ from the absolute central frequency of each line. The 
frequency spacing used here is the same as the spacing of data points '{DFLTAG 
in the DERPTH program) in the Venus A and Venus 5 planetary spectra, and 
the total width of the region over which the absorption is calculated is 
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large enough to include the most significant parts of the profiles of the 
HCJt lines f whose half-widths in the planetary spectra lie between 0.01 
and 0.02 cm FORTRAN routines written by Armstrong '(1967) were used to 
evaluate numerically the integral in Equation (4.10) . Then, after the 
calculation of the absorption coefficient for all levels. Equation (4,18) 
was integrated numerically using the trapezoidal rule to obtain the trans- 
mittance at each frequency for the atmospheric layer between 95 km and each 
level of interest. Since the Venus atmosphere is very thin above 95 3on, 
with total pressures less than 0.01 mb, the contribution to the transmittance 
arising' from layers above 95 km is negligible when compared with the con- 
tributions from lower layers. Thus the calculated values for the layer up 
to 95 km can be considered to represent the total transmittance from the 
level of interest to the top of the a'tmosphere. As the final step in the 
calculation of realistic HCS, line profiles, the theoretical transmittances 
were convolved with the ILS for the unapodized spectra. 

Figures 23 and 24 show a sample of the results for the median model, 
giving profiles of the absorption coefficient (Pigtire 23) and transmittance 
(Figure 24) in the region of the R4 line. In the four graphs the 

letters plotted refer to different significant levels in the model atmo- 
sphere, which are also indicated on the vertical axis in Figure 22. The 
absorption coefficient profiles for the different levels shown in Figure 23 
indicate that near the peak of the line, the maximum contribution to the 
absorption comes from the upper layers of the atmosphere, while near the 
base of the line the maximum contribution is from lower atmospheric levels - 
Figure 24 illustrates the net transmittance between the top of the atmo- 
sphere and each of the significant levels. These significant levels can be 
viewed also as the locations of the reflecting surface at the top of the 
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cloud deck. The most realistic locations for this surface are between 60 
■and 69 km;. absorption and transmittance for these layers are represented in 
the figures by curves B through E. Note that there is increasing asymmetry 
in the line shapes for the lower > regions of the atmosphere {curves E, F, 
and G) . This is due to the increasing negative displacement of the cen- 
tral peak of the absorption coefficient profile with depth into the atmo- 
sphere. Careful comparison of the line shapes calculated separately using 
the Venus A (Figure 24a) and Venus 5 (Figure 24b) velocity scaling factors 
reveals subtle differences in the two as 3 ^mmetric line shapes for the same 
model atmosphere. This is because the identical profile of retrograde zonal 
velocities produces a positive Doppler shift for Venus A and -a negative one 
for Venus 5; coupled with the pressure-induced lineshift, which is always 
negative for the R4 line, the difference in sign of the' Doppler shifts causes 
the Venus 5 line profiles to be even more skewed toward smaller frequencies 
than would be expected for the Venus A profile. 

In Figure 25 there is reproduced the R4 line profile for trans- 

mittance between 63 and 95 km in the Venus 5 median model atmosphere, which 
is the same as curve "D" in Figure 24b. The five equally-spaced chords 
were drawn across the line profile and their central positions calculated 
using the same algorithms as those in the DERPTH program- The heights in 
-the model atmosphere to which these chord-center shifts correspond are given 
along the right-hand vertical scale of the figure. The details of calcula- 
tion of these heights from the shifts are given in Section 4.4.1. The 
specific distribution of the chord-center shifts with height in a spectral 
line depends on the particular line, as well as on the pressure and wind 
velocity profiles. For the line shown in the figure, the shifts are all 
negative, and their magnitudes increase monotonically from the peak to the 


base of the line. 
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The question arises whether. the correspondence between each equally 
spaced chord and a given atmospheric pressure level remains the same for all 
lines in the absorption band. . Shifts of the five chord centers as shown 
in Figure 25 were calculated for each of the eight lines in the HCJl R-branch 
and the heights/ pressures and wind velocities corresponding to these 
chords were compared. Results are given in Table 6 for the A-H{1) model 
velocity profile for both Venus A and Venus 5 synthetic line profiles. In 
the table are listed the heights in km in the model Venus atmosphere cor- 
responding to each chord for the individual lines/ along with mean values 
and standard deviations for the whole R-branch. 

There appears to be more variation in the reference heights for 
chords Cl and Cz relative to the other chords. This variability is due to 
the fact that the spacing of data points in the line profile is 0.002 cm 
allowing calculation of the chord-center positions within only ±0.002 cm 
For lower J values and near the peaks of the lines/ the magnitudes of the 
pressure- and wind- induced lineshifts are generally less than about 0.001 
cm so that the uncertainty in determining the positions of the chord 
centers is very significant. For chord's C 3 and c^, in the middle and lower 
portion of the absorption line, most of the expected lineshifts are consider 
ably larger (from 0.001 to about 0.005 cm ^) , and the uncertainty is not so 
important. The reference heights for chord cs also show considerable vari- 
ability due to noise in the line profile near the mean background level (FM 
in Figure 25} caused by the side lobes of the instriimental line shape for 
the unapodized spectra. The maximum- amplitude of these side lobes is about 
0.18, so that the noise introduced into the line profile by the side lobes 
will extend to a level of at most 0.82, which is less than the transmittance 
at the level of chord cs for all blit the strongest lines in the band. 

Also in Table 6 it appears that the reference heights for the chords 
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in lines having lower J numbers are somewhat greater than those determined 
for the other lines in the band. This effect is an artifact of the method 
by which the chords are drawn across the line profile. The DEBPTH program 
studies the symmetry of each line over a constant interval of width 6*DS, 
where DS is a programmer- specified average half-width for all the lines of 
interest in the spectrum. As illustrated in Figure 17, the program finds 
the difference between the mean backgroxmd and minimum intensities within 
this interval (DGMAX = FM - GGMIN) , and takes 1/6 of this difference as the 
spacing for the chords. If the interval 6*DS is narrower than the width of 
the absorption line at its base, the mean background level will be drawn at 
a higher level in the line profile, and the five chords will be drawn only 
in the upper part of the line. In the model calculations given in Table 6, 
the width of the spectral interval 6*DS was 0.200 cm which caused the 
background level FM to be located at heights less than 0.99 in the "D“ line 
profiles for the RO through r 3 lines. A survey of the actual halfwidth of 
the HCA lines in the planetary spectra shows that for Venus A the largest 
halfwidth is 0.0138 cm while the interval 6*DS is 0.1020 cm and for 
Venus 5 the largest halfwidth is 0.205 cm ^ and 6*DS = 0.3660 cm Thus 
we expect all the hcA line profiles in the planetary spectra to be com- 
pletely contained in the interval 6*DS, and that the mean reference levels 
for the low-J lines calculated from these spectra will not differ signifi- 
cantly from the mean values for the whole band. Since the half -widths for 
the calculated "D" line profiles are close to those measured for the HCfi, 
lines in the spectra, the standard deviations given for the reference levels 
in Table 6 are probably not very different from those for the planetary spectra. 
Barring other sources of error, we would expect to be able to calculate ref- 
erence heights within at least ±1.5 km, which corresponds to an error of at - 
most 30 percent in atmospheric pressures in the region 60 to 70 km above the 
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planet's surface. 

The model results for shifts of the central positions of the chords 
for all the HCA lines are presented graphically in Figures 26 through 28. 

In both the A-H(l) and median models, the velocity is increasing with depth 
into the atmosphere. Thus not only do the corresponding lineshift curves 
have increasing slopes, reflecting pressure- induced lineshifts increasing 
with depth, but also the entire curve for each level is displaced due to 
the wind-related Doppler shift. This Doppler shift can easily be found from 
the graphs, by measuring the displacement of the zero-crossing point for 
each curve, which is located between the R2 and R3 lines at approximately 
5730 cm Note that for the A-H(2) model (Figure 28), which has a con- 
stant velocity at all levels, the displacements of the zero-crossing points 
are identical, and all the pressure-shift curves pass through the same point. 
Although this may also appear to be the case in Figures 26 and 27, close 
examination of the shifrs shows that the curves do not intersect at exactly 
the same point for these two cases where the wind velocity is not constant 
with height in the atmosphere. 

The zonal velocity profiles retrieved from the calculated line pro- 
files in all four models are shown in Figure 29 along with the original wind 
profiles for each model. With the exception of model A-H(2B), which has 
sharp variations in zonal wind velocities within a few kilometers, the 
retrieved profiles appear to closely approximate the original ones. A more 
detailed discussion of the results and sources of error in these retrieved 
model profiles, along with the results for the planetary spectra, will be 


undertaken in Section 4.4.2. 
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4,4 The inverse problem; estimation of 
the velocity profile from measured Ixne ' 
asymmetries in the HCiZ. 2~0 band 

In the previous section it was shown by direct calculation how 
pressure-induced lineshifts and wind-induced Doppler shifts can be respon- 
sible for detectable asymmetries in ECZ spectral lines in high-resolution 
planetary spectra. It was also shown that the shifts in the central posi- 
tions of chords equally spaced from the base to the peak of each line con- 
sistently represent the effects of wind velocity and pressure at distinct 
atmospheric levels for the RO through R6 lines in the HC5. 2-0 band. Thus 
the use of measured line asymmetries to calculate both atmospheric pressures 
and zonal wind velocities at specific levels seems to be a viable means of 
estimating vertical profiles of wind velocity in the Venus atmosphere above 
the clouds. 

4.4.1 Procedure 

The' DERPTH program calculates and prints out for each spectral line 
the central position determined for chord ci (Figures 17 and 25) and the 
deviations of the central positions of the four succeeding chords from that 
of Cl. This program was used for the Venus 5 averaged spectirum and for all 
of the individual spectra. An earlier version of the same program called 
DEEPTX was used for the Venus A averaged spectrum; in this case nine chords 
were drawn across each spectral line and had their central positions determined. 

In each spectrum all the lines in the R-branch of the and E^^CZ 

2-0 bands having quality factors of 3 or greater were used in the calculation 
of the pressures and zonal wind velocities. For each line there first was 
calculated a Doppler-corrected absolute position/ i.e./ 

Av 


= (1 + (■ 


(4.19) 
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where u is the absolute vacuum wavenumber for the line measured by Guelachvili 
o 

(1976) and Av/v is the combined translational and earth-rotational Doppler 
shift factor calculated from Ephemeris data for the spectrum. The Doppler- 
correct^ absolute position was then subtracted from the central position 
determined by the program for chord ci (listed in Table 4) in order to ob- 
tain the lineshift for that chord, 5i. Lineshifts for the succeeding chords 
were then found by adding the value of 5i to the chord-center displacements 
(62 ” 5x, 63 - 5i, etc.) which had been previously calculated by the DERPTH 
or DERPTX programs. The resulting lineshifts for all the chords are plotted 
for the whole R-branch in Figure 30 for the two averaged spectra Venus A and 
Venus 5. The points plotted in the figure are averages of the shifts for 
the and lines, weighted by their respective guality factors; 

this averaging produces a "good" data point for each J value. Chord-center 
shifts were also measured in individual spectra recorded on the first and 
last days of observations included in these averages. The patterns of the 
lineshifts for the individual spectra are quite similar to those for the 
.averaged spectra. 

An initial value of the wind-induced Doppler shift component 
of the lineshift for each chord was estimated from the graphs. Pressure- 
shift factors r^ were then calculated for each line in the R-branch, using 
the formula 


r . 
1 



(4.20) 


where 6^ is the net shift of the chord center, and 5^ is the laboratory 
measured lineshift averaged for both isotopic species, which refers to a 
standard CO 2 pressure of 122.2 mb. 
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Several iterations of this calculation of r. for the R~branch lines 

1 

were necessary in order to find the best value of 'tfie estimated 

value for is too large and negative, the calculated values of r^ will 
be small for the higher-J lines relative to those for the RO and Rl lines; 
similarly, a too-large positive estimate for will cause the RO and Rl 
lines to have excessively high values of r^. The iteration of the estima- 
tion of and the calculation of r, continues until the values of r, ob- 
Di 1 1 

tained • appear to be evenly distributed with respect to J. The final values 
of r. are then averaged over J and multiplied by 122.2 to obtain the average 
pressure at level i in mb- 
it is assumed that the wind-induced Doppler shift 6 ... refers to the 

Di 

displacement of the zero-crossing point of each pressure-shift curve, which 
occurs at about 5730 cm ^ . Thus v . , the spectrally-determined earthward 
component of Venus zonal velocity can be found using the formula 


V , 

si 



5730 cm 


* c 


(4.21) 


where c is the velocity of light. The equatorial zonal velocity can then 
be found through dividing v^^ by the velocity scaling factor derived in 
section 4.2.2. 

4.4.2 Results and error analysis 

Velocity profiles retrieved from the calculated model line pro- 

files are shown in Figiire 29, and the results obtained from the planetary 
spectra are given in Figure 31. The vertical error bars on the points plotted 
in Figure 29 represent the standard deviation from the mean reference pres- 
sure determined for each chord, averaged over the RO through R8 lines of the 
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2r0 band. .-The standard- deviation- is- about iS percent for chord ci , 

J.5 .percent for chords C 2 ~and-C 3 ,Tandr-20* toi_2S percent for chords ct, and C 5 . 
This result is. consistent. with thei variation of-the-reference heights for the 
chords! in. model Ar-H(l) presented. in-Table ~ 6 and discussed in Section 4.3. 

Inr the.' results for-.the planetary, spectra the standard deviations in the cal- 
culated pressure levels, are- approximately, twice-' as-lafge- as those for the 
model; line profiles. --This -additional uncertainty^ is most probably due to 
variations in the„HCJl line profiles 'in the planetary" spectra resulting from 
overlapping of those lines with- other. absorption features arising in the 
atmospheres o£ both Venus and the earth. Such super imposition of absorption 
features. was not accounted for when- the -model- line profiles were calculated. 

: . ■ . From the model calculations we can also estimate a minimum error 
limit on the zonal wind velocity determined for "each pressure level. The 
displacement of the zero-crossing point of the pattern of lineshifts 

for each chord can be estimated graphically within ±0.05 x 10 cm , so that 
at best the corresponding zonal wind velocity can be determined within ±3 m 
sec Since the displacement of the zero-crossing point affects the esti- 
mate of pressure for each chord-center shift, there exists a direct relation- 
ship between the uncertainty in the reference pressure and uncertainty in 
the wind velocity calculated for each chord. In both the model calculations 
and the planetary spectra results, it appears that for chords near the peaks 
of the HCA lines, a 1 -percent error in estimation of the pressure corresponds 
to about a 2-percent error in the zonal wind velocity at that level. The 
range -of possible wind velocities corresponding to a given pressure seems 
to increase with depth into the atmosphere; for chord C5, near the base of 
the line and heavily influenced by background noise, a 1 -percent error in 
calculated pressure can correspond to as much as 25 percent error in the 
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calculated wind velocity. In other words, for the profiles derernu-ned from 
the planetary spectra, the uncertainty in the estimated zonal wind velocity 
ranges from about 30 percent in the upper part of the profile to several 
hundred percent for some values estimated from chord C 5 . Error bars showing 
the variability of the calculated wind velocities are shown in Figure 31; 
in the interest of clarity they are drawn only for the two profiles obtained 
from the averaged spectra. 

One additional comment must be made concerning the model results. 

In model A-H(2B) the original velocity profile showed an abrupt decrease in 
wind speed from 130 m sec ^ to zero in a narrow layer about 2 km thick. How- 
ever, the retrxeved velocity profile shows a much more gradual change in 
wind velocity with height than was actually the case. This suggests that 
retrieval of wind velocity profiles from pressure-induced lineshifts is rela- 
tively insensitive to variations in the true velocity profile with vertical 
scales ' smaller than about 5 km. However, as can be seen from the results 
for the other three models presented in Figure 29, the mean shape of the ve- 
locity profile can be detected quite well"; It may be a worthwhile future 
endeavor to calculate the shape of weighting functions for the retrieval of 
wind velocxties from measured asymmetrxes in various parts of a spectral line 
profile - 

Velocity profiles estimated from the two averaged planetary spectra 
(solid lines) and from four of the individual spectra (dashed lines) are 
shown in Figure 31. It is apparent in the figure that the profiles obtained 
using rhe individual spectra have generally the same shape as those obtained 
from the corresponding averaged spectra, but the magnitudes of the velocities 
vary greatly from one profile to the next. It was noted in the preceding 
section that the patterns of the chord shifts are similar for each averaged 
spectrum and its component spectra. Generally the wind velociries determined 
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from the individual spectra appear to fall within the range of velocities 
given by the error limits for the averaged spectra. However, careful com- 
parison of the individual shift patterns also shows that the wind-induced 
Doppler displacements of the pressure-shift curves are generally more posi- 
tive on the first day of observations and more negative on the last day. 

-This suggests that there may be systematic errors in the calculation of the 
translational and/or earth-rotational Doppler shift factors, probably as- 
sociated with the averaging of direction and distance of Venus from the ob- 
servation site diiring the 2 to 3 hours necessary for each recording. An 
alternative explanation for these apparently uniform differences between 
velocity profiles from individual spectra is that they may actually represent 
real changes in the zonal velocity field during each period of observations 
(10 days for Venus A, 5 'days- for Venus 5) . 

An independent estimate of zonal wind velocities in the upper atmo- 
sphere of Venus was made from the same set of spectra, using lines in the 
qIZoIS 2 _g band near 2.35 ym. Nearly all the lines in this band from P24 
through R24 appeared in the planetary spectra, and a little more than half 
of these were sufficiently unblended with other absorption features to be 
useful in calculating the Doppler shift factor Av/V for the spectrum. In 
these calculations the absolute wavenumber measurements of Guelachvili (1973) 
were used as standards for the spectra in the Venus A series, while for the 
Venus 5 series the measured positions of lines arising from the CO cell which 
had been placed in the path of the interferometer were used whenever possible. 
According to the results of L. Young (1972) and Connes et al . (1968) , the 
mean pressure of line formation for CO in the Venus atmosphere is 36 to 60 
mb, corresponding to a height of 67 to 70 hm above the surface. At these 
pressures the magnitudes of lineshifts induced in CO, by CO 2 are likely to 
be no greater than 5 x 10 cm ^ , 


and can be neglected. The difference 



58 


between the Doppler shift factor calculated from the CO 2-0 band and that 
calculated from the Ephemeris data for each spectrum represents the Venus 
wind velocity at the mean level of CO line formation. The results are 
listed in Table 7. It can be seen that the velocities calculated from the 
CO band generally are consistent with those calculated from the upper parts 
of the HCil lines in the same spectra- Nevertheless there again appear the 
wide variations in velocities determined from different individual spectra 
in the same series, which adds to the evidence that these variations are 
caused either by averaging errors in the calculation of the Doppler shift 
factors from the Ephemeris data, or by actual changes in the zonal velocity 
field. 

For comparison to the velocities estimated here from HCil lines in 
the infrared spectra, the maximum zonal winds measured by Veneras 4, B, 9 
and 10 between 40 and 60 km are also indicated in Figure 31, The shaded 
region in the figure represents the range of velocities estimated from UV 
photographs and other spectroscopic studies. It is noted that both of the 
velocity profiles estimated from HCA lines in the averaged spectra contain 
some values that lie in this shaded region, which indicates that the magni- 
tudes and directions of these estimated velocities are more or less realistic. 

Another estimate of the range of realistic wind velocities possible 

in the Venus upper atmosphere can be made by calculating the speed of sound 

waves for this region, using the relation c (sound) = , Zonal wind 

CO2 

velocities would not be expected to exceed this magnitude. For the tempera- 
tures in the Venus atmosphere between 60 and 80 km, the calculated speed of 
sound is about 250 to 300 m sec Thus, the velocities estimated from the 
two averaged spectra appear to be well within realistic limits. 

The velocity profile estimated from the Venus A spectrum shows a wide 
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variation in zonal velocities within an atmospheric layer only 10 3on in depth. 
Such shears may not be unreasonable, however, since as Limaye (1977) has 
podLnted out, the atmosphere is stably stratified in this region, and strong 
vertical mixing is unlikely. Also, the scale height in this region is 
relatively constant at 4 to 5 km. Near 60 km there appear to exist strong 
retrograde winds of about -75 m sec ^ , and there is slight evidence for a 
relative maximum in the profile of velocity magnitudes at that height. The 
retrograde winds decrease rapidly with height, becoming zero at roughly 64 
km. Above this height the winds are apparently prograde in direction, with 
magnitudes steadily increasing with height to a maximum of about +68 m sec ^ 
near 70 km, where the estimated profile ends. The detection of prograde 
wind velocities with magnitudes of several tens of meters per second is in- 
deed a remarkable result. The only other measurements of prograde winds are 
several single spectroscopic estimates by Traub and Carleton (1975) where 
the magnitudes of the prograde winds were only a few meters per second. 
Nevertheless, as noted in the discussion in Chapter III, some plausible dy- 
namical regimes for the Venus atmospheric circulation, notably the model of 
Ramanathan and Cess (1975) , do admit the possibility of prograde winds in 
the upper part of the Venus stratosphere. 

The Venus 5 velocity profile appears to refer to a shallower layer 
in the atmosphere. In this region (about 65 to 70' km) all the estimated 
velocities are retrograde, having magnitudes from 90 to 131 m sec The 
Venus 5 velocities appear to increase in magnitude with height, showing con- 
siderably less vertical wind shear than in the Venus A profile. 

Since the Venus A ’and Venus 5 velocity profiles, as well as the other 
estimates of wind velocities in the Venus atmosphere up to 80 km, refer to 
different local times of day, it is helpful to draw a cross section in 
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, longitude and -height of all available measured Venus equatorial zonal 
_velocities. . Figure 32 is such- a cross section; the -plotted data points 
^include velocities from the profiles 'estimated front the Venus A -and Venus 5 
^ispectra, along with measurements from the Venera' probes/ and' the veloci- 

estimated by Limaye (1977) and Traub and Carleton -(1975)-^- Each veloc- 
ity-measurement is .piotted at approximately the mean height and longitude 
(relative to. .the. sub-solar point) to which :if refers; error'-'bars are not 
■included for the sake of a clearer drawing. Velocity magnitudes are iso- 
.plethed using dashed lines, and the solid lines with arrows indicate pos- 
sible streamlines for the mean flow. . . - -- - - - 

It is important to remember that the zonal velocities plotted in 
the figure do not represent simultaneous observations, or, except for the 
Traub and Carleton' results, averages over time periods longer than about 
- a week. Thus this cross section may not necessarily represent the steady- 
estate 'equatorial circulation of. the Venus' atmosphere. According-- to Dollfus 
(1975) , in his review of photographs of the planet taken between 1927 and 
.1974, the general circulation of the Venus atmosphere as represented by 
'r the motions of large UV features has" remained more or less- the same within 
that. 47-year time period. However, one cannot rule out the possibility of 
changes in the circulation occurring below the level of the UV-reflecting 
layer which may have time scales of the order of months or years. One 
short-term phenomenon detectable from Earth rs the brighrening of the polar 
regions which occurs on a time scale of several months and may be a "sea- 
sonal" effect. As Dollfus (1975) has noted, the brightening events occur 
almost exclusively when the planet is in the same half of its orbit around 
the sun. It is not yet known wherher. this brightening of the polar caps 
is related to some subtle ^ change in the circulation of the Venus atmosphere. 



61 


Long-term changes in the Venus atmospheric circulation may be re- 
sponsible for the inconsistency of the flow field due to the Venera 4 and 
Venera 8 velocity profiles in relation to the' circulation derived from 
the other measurements. These two Venera profiles were measured prior to 
1970, while the other results plotted in Figure 32 came from the period 
from 1971 to 1976. The large mass transports associated with the high 
zonal wind speeds near 40 Ion at the morning terminator reguirsv. either ex- 
cessively high zonal wind speeds (>1000 m sec in the stratosphere or 
meridional velocities several orders of magnitude larger than have been ^ 
obsejrved in order, for the mass balance in the general circulation to be 
maintained. Because of this inconsistency and the lack of additional ob- 
servations, we shall attempt no further discussion of the tropospheric cir- 
culation. except to say that the zonal winds in the Venus troposphere appear 
to be consistently retrograde in direction. 

Given the sparseness and relative uncertainty of the data points, 
as well as the weak mass transport in the stratosphere, some features of 
the circulation pattern can be resolved. In the layer roughly from 65 to 
70 km, prograde velocities in the afternoon longitudes and retrograde winds 
in the morning region result in a zone of convergence at the sub-solar point 
and slightly westward, producing net upward vertical .motions above this 
layer and downward motions .below 'it . Above 75 km the flow appears to be 
uniformly retrograde in -direction. 

The stratospheric velocity field shown- in Figure 32 is suggestive 
of a two-celled bipolar circulation schematically similar to that generated 
by the "moving flame” mechanism. However, more complex diurnal circulation 
patterns may show the same type of flow field in the sub-solar region. For 
example, in the three-dimensional model of Young and Pollack (1977) , where 
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the stratospheric circulation is generated by the nonlinear interaction of 
planetary-scale eddies with the mean Hadley circulation of the troposphere, 
there also results net upward motions in the s'ub-solar region and consis- 
tently retrograde zonal winds at higher altitudes. It is not possible to 
determine from the sparse data in Figure 32 exactly what type of mechanism 
may be responsible for the stratospheric circulation. 

Additional evidence for strong vertical motions in the sub-solar 
region was found in the Mariner 10 ultraviolet photographs of the equa- 
torial regions. Belton ^ al . (1976b) have identified a region called the 
sub-solar disturbance (SSD) which contains cellular features having linear 
scales of 200 to 500 1cm, indicating the presence of convective motions in 
this region. The cellular features do not appear upwind or far downwind 
of the sub-solar point, where the UV features take the form of east-west 
streaks indicating retrograde zonal motion - 

The retrograde velocities above 72 km in Figure 32 correspond with 
the pattern of UV features found in Mariner 10 photographs. However, none 
of the photographs shows obvious evidence for prograde zonal winds. The 
prograde velocities occurring at 65 to 70 km in the figure would lie be- 
low the height estimated by Limaye (1977) for the UV-reflecting layer and 
thus would not necessarily be detected in DV photographs _ 

In summary, it appears that the velocity profiles determined from 
HC£ lines in the Venus A and- Venus 5 spectra appear to be consistent with 
other estimates of zonal wind velocities and with the patteras of UV fea- 
tures photographed by Mariner 10. The composite stratospheric equatorial 
wind field appears to represent a diurnal ly-yarying circulation induced by 
local heating in the sub-solar region. The exact mechanism responsible for 
this circulation cannot be determined from the currently available data. 
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CONCLUSION 

The results presented in this study allow us to draw" a number of 
conclusions; some of these have been stated in earlier chapters, but will 
be restated and summarized here, 

(1) Pressure-induced lineshifts for ,HCA and HF absorption lines due to 
the CO 2 iu the Venus atmosphere can be detected and measured in high- 
resolution infrared spectra of the planet. The lineshifts should be de- 
tectable in any infrared spectrum where the line positions can be determined 
within ±0.001 cm Since spectra obtained in narrow frequency regions 
using tunable diode lasers, as well as spectra from Fourier- transform inter- 
ferometers, can achieve such accuracies, the possibility exists for develop- 
ing suitable monitoring instruments making use of pressure-induced line- 
shifts for hydrogen halides perturbed by CO 2 (for the Venus atmosphere) or 
by N 2 and O 2 for the earth's atmosphere) to measure wind velocities and 
pressures in planetary atmospheres. Since the lineshifts are linearly de- 
pendent on atmospheric pressure, they can be used to determine the pressure 
of line formation for the hydrogen halide lines in the planetary atmosphere, 
and Doppler shifts due to atmospheric motions at the pressure levels deter- 
mined from the lineshifts can also be resolved .and .measured. 

(2) For the Venus spectra used in this work, it was shown that asym- 
metries in line shapes in the HC£ 2-0 band, caused by both pressure-induced 
shifts and Doppler shifts due to atmospheric motions at different levels in 
the upper atmosphere, can be used to estimate vertical profiles of the 
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average equatorial zonal wind velocity as a function of Venus atmospheric 
pressure. The R-branch of the HCA 2~0 band in the Venus spectra was chosen 
mainly because of the excellent quality of the spectral lines in this region. 
However, vertical profiles of velocities in planetary atmospheres can also 
be estimated using asymmetries of absorption lines belonging to other mole- 
cules, provided that atmospheric spectra of sufficient quality exist, and 
that there are available' complete laboratory measurements of the pressure- 
induced iineshifts of interest. 

' - In' the case of' Venus, the atmospheric pressures determined from the 

HCi line asymmetries are reasonably accurate, but the corresponding wind 
velocities show considerable uncertainty. A major portion of the uncertainty 
may be due to' error in calculation of the translational and earth-rotational 
Doppler shifts. This" error arises from the fact that the interferograms 
were not meastared simultaneously; each one was integrated over a two- to 
three-hour period*, and the final spectra, Venus A and Venus 5, are averages 
of several recordings. 

An ideal experiment, focusing on the accurate measurement of line 
asymmetries in order to 'obtain Venus wind profiles, would involve the fol- 
lowing: (i) very high- resolution measurements in the region of the ECS. 2-0 

band only, with as short a sampling time as possible; (ii) measurements made 
in pairs, one with the viewing aperture centered at the sub-earth point, where 
the earthward component of zonal velocity is zero, and the other made at ap- 
proximately the same time with the viewing aperture centered on the point 
of interest on the planetary disk- This scheme is similar to that used by 
Traub and Carleton (1975) in their spectroscopic measurements of wind veloc- 
ities on Venus. The instantaneous combined translational and earth-rotational 
Doppler shift could then be determined directly from the reference spectrum. 



65 


eliminating much of the averaging errors encountered in the present work. 

(3) Venus zonal velocity profiles estimated from the infrared spectra 
show large- vertical wind shears in the region 60 to 70 kra above the surface 
of the planet, and the velocity profiles appear to vary greatly with loca~ 
tion relative to the sub-solar point. The profiles obtained from the Venus 
A spectrvmi, which refers mainly to local afternoon regions on 'the planet, 
show a dramatic change from retrograde winds near 60 km to prograde veloc- 
ities of the order of several tens of meters per second near 70 km. When 
these profiles are considered along with the few existing measurements of 
Venus wind velocities above 50 km, the composite stratospheric wind field 
appears to represent a diurnal circulation induced by strong local heating 
in the sub-solar region. Such a circulation pattern also appears to be 
consistent with the configuration of UV features photographed by Mariner 
10 in the equatorial regions. A definite identification of the mechanism (s) 
responsible for the observed stratospheric circulation cannot be made* at 
this time due to the uncertainties in existing measurements and the absence 
of wind, data for the night side of Venus. 

The results of this study, besides demonstrating how- asymmetries 
in absorption lines can be used to estimate zonal velocity profiles in the 
Venus atmosphere, have provided some additional information regarding the 
vertical and longitudinal variation of zonal winds in the Venus stratosphere. 
-It is hoped that Doppler tracking of the descent probes of the Pioneer 
Venus Mission in December 1978 will provide much more detailed information 
on the velocity distribution in the Venus atmosphere. The four probes are 
preprogrammed to impact the planet at four widely separated sites, includ- 

< 

ing two on the night side, and are expected to provide a nearly simultaneous 
sample of Venus wind velocities between the surface and 70 km. Perhaps the 
additional information will verify the results obtained in the present work. 
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TABLE 1 

PHYSICAL CHAEACTERISTICS OP VENUS AND E2iRTH 

\ 




Venus 

Earth 

Equatorial radius (km) 




(Radar) 

* 

6052 ± 5 

6378.2 

(Optical) 


6115 ± 13 


Mean distance to Sun 




(million km) 


108.6 

149.5 

(a.uO 


0.723332 

1-000000 

Mass relative to that 
of Earth 


0.81427 

1 

Acceleration due to' 
gravity (cm sec ^) 


885.15 

980.1 

Solar constant (ly min ^) 


3.8 

2.0 

Bond albedo 


0.77 ± .07 

0.3 

Inclination of spin axis 
to. orbital plane 


3 

23.5 

Inclination of orbit to 
the ecliptic 


3.394 

0 

Eccentricity of orbit 

' 

0.0068 

0.0167 

Length of year (days) 


225 

365.26 


Period of rotation of 

solar planet (days) 243 (retrograde) 1 

Tangential speed of solid 
planet due to rotation 

(m sec ^) 1.8 cos (j) 463.8 cos 

Coriolis parameter 
(10 ^ sec ^) 


5.99 sin ({) 


1454 sin (j> 



TABLE 2 


COMPOSITION OF THE VENUS ATMOSPHERE 




(a) Molecules observed in the 

atmosphere of Venus 

above the clouds 

G3.S 

Mixing Ratio 

1 Total Amount 

(cm-atm___ ) 

' STP 

Pressure 

(mb) 

Temperature 

(K) 

Reference 

CO 

4.5 ± 1.0 X 

lO”® 13 

60 

240 

Connes et al. (1968) 


5.1 ± 0.1 X 

lo"® 20 

36 

249 ± 3 

L. Young (1972) 

CO 2 

0.97 ± 0.04 

- 

- 

- 

Vinogradov et al. (1970a,b) 


- 

3.3 ± 0.3 X 10® 

100 

240 ± 10 

Connes ^ al. (1967) 


- 

2.8 ± 0.4 X 10® 

60 

240 

Connes ^ al. (1968) 


- 

4.1 ± 0.1 X 10® 

36 

249 ± 3 

L. Young (1972) 


- 

0.9 - 5.5 X 10® 

100 

245 

Fink et al. (1972) 

mi 

6 X 10 ’ 

1.9 ± 0.4 X 10~^ 

80 

270 ± 30 

Connes et al. (1967) 


4.2 ± 0.7 X 

1.7 ± 0.3 X lo"^ 

160 

240 ± 16 

L. Young (1972) 

HF 

5 X 10“® 

- 

100 

240 

Connes et al. (1967) 


7 X lo"® 


67 

249 ± 3 

L. Young (1972) 



10““ 

35 

2.5 X 10““ 

100 

8.7 X 10“® 

25 

io““ 

39 

io"“ 

30 

1.1 X 10 ^ 

.. 

2 X 10“® 

0.6 

2 X 10“® 

5 

10“ ® 

t 

0.6 X 10“® 

0.2 


0.5 


90 

600 


600 


10000 


- Dollfus ( 1963a, b) 

Bottema ^ (1964a, b) 

^ M 

- Belton & Hunten (1966) 

Spinrad & Shawl (1966) 

295 Vinogradov ^ (1968) 

Kuiper (1969) 

Schorn ^ al. (1969a, b) 

500 Vinogradov et al. (1970a, b) 

-- Pink et.al, (1972) 


Gull et al. (1974) 



TABLE 2 — Continued 


(b) Molecules not observed in the atmosphere of Venus above the clouds 


Gas 

Maximum 

Maximum amount 

References 


Mixing Ratio 

(cra-atm„„„) 

STB 


CH 2 O 

10“ ® 

I 

0 

Owen and Sagan (1972) 

CH 3 and 
aldehydes 

10“ ® 

10"^ 

Owen and Sagan (1972) 

CH 4 

10"® 

3 X 10"^ 

Connes et al. (1967) 

CHaCil 

10“ ® 

3 X lO**^ 

Connes et al. (1967) 

CH 3 F 

10"® 

3 X 10"^ 

Connes ^ al. (1967) 

CH 3 COCH 3 and 
ketones 

!-• 

0 

1 

o> 

10”^ 

Owen and Sagan (1972) 

C 2 H 2 

10 ~® 

3 X 10“^ 

Connes ^ al. (1967) 

C 2 H 4 

10 ~® 

2 

Kuiper (1952) 

C 2 H 6 

10”® 

4 

Kuiper (1952) 

C 3 O 2 

10“^ 

10“^ 

Owen and Sagan (1972) 



HCN 

10 "® 

1 

3 X 10"‘ 

'1 1 1 i 

Connes ^ al. (1967) 

HaS 

10 "^ 

10 “^ 

Owen and Sagan (1972) 

Na 

2 X 10 "^ 

-r' 

Vindgradov ^ al. (1970a, b) 

NHs 

3 X 10 "® ' 

10~^ 

1 ‘ 

Kuiper (1962) 

NO 

10~® 


Owen and Sagan (1972) 

NOa 

io“® 

10 "® 

Owen and Sagan (1972) 

NaO 

6 X 10“** 

200 

Kuiper (1952) • ■ t > 

Na04 

4 X 10“® 

4 X lo"® 

Owen and Sagan (1972) 

Oa 

7 X 10 "® 

24 

Belton and Hunten (1968) 

O 3 

3 X 10“® 

3 X 10 “'* 

Owen and Sagan (1972) 

SOa 

5 X 10 "® 

5 X 10 "® 

Cruikshank and Kuiper (1967) 


lo“® lo"^ 


Owen and Sagan (1972) 



TABLE 3 


PARAMETERS FOR RECORDED SPECTRA OP VENUS 



Spectrum Date Times (Local) 


Illuminated Disk Distance from Semidiameter 

Earth (a.u.) (arc sec) 







Start 

End 

Phase 

Angle 

Fraction 

Ilium. 

Stellar 

Mag. 




143 


14 

Oct 

73 

14:28 

17:30 

74.9° 

0,631 

-3.8 

0.897 

3169 

9.37 

149 


15 

Oct 

73 

14:12 

16:26 

75.3° 

0.627 

-3.8 

0.889 

8692 

9.45 

156 


22 

Oct 

73 

15:13 

17:45 

78.6° 

0.599 

-3,8 

0.837 

5038 

10.04 

161 


25 

Oct 

73 

14:49 

17:50 

80.0° 

0.586 

-3.9 

0.814 

9448 

10.32 

Venus 

A 

average 




77.2° 

0.611 

-3.8 

0.859 

9087 

9.80 

226 


12 

Jun 

74 

4:48 

7:27 

57.4° 

0.769 

-3.5 

1.201 

0340 

7.00 

228 


13 

Jun 

74 

5:10 

7:41 

57.0° 

0.772 

-3.5 

1.207 

8121 

6.96 

230 


14 

Jun 

74 

4:29 

7:37 

56.6° 

0.775 

-3.5 

1-214 

5570 

6.92 

231 


15 

Jun 

74 

4:29 

7:40 

56.2° 

0.778 

-3.4 

1.221 

2683 

6.89 

234 


16 

Jun 

74 

4:30 

7:38 

55.8° 

0.781 

_3.4 

1.227 

9454 

6.85 

Venus 

5 

average 




56.6° 

0.775 

-3.5 

1.214 

5234 

6.92 


-a 

H* 


TABLE 4 


HYDROGEN HALIDE LINES OBSERVED IN THE VENUS SPECTRA 


(a) HF 1-0 band 


Line 

Standard position (cm 

1 

Venus A 

* ) ' i 

i • 1 

Venus 5 

‘111 ' 111 1 • 1 1 1 li I'*' { 

• w 


(±2.5 X 10 ** cm 

Central Position 
(cm“M 

1 1 . 1 1 ‘ < 1 

Quality' 

Factor 

J 

•' Central Positioh ' 

(cm 

Quality 

Factor 

R1 

4038.9625 

4039.135106 

1 

4038.807469 

3 

R2 

4075.2936 < 

'4075.467221 ' ‘ 

0 

" '' 4675.138583 “ ' 

1 

R3 

4109.9363 

4110.111393 

3 

' 4109.776186 ’ ‘ 

1 

R4 

4142.8460 

4143.020316 

1 

’ 4142.685340 *" 

2 

R5 

4173.9796 

4174.155088 

3 

4173.818178 

1 

R6 

4203.2960 

4203.475137 

0 

4203.1315 

0 



(b) HP 2-0 band' 





(±15 X lo”** cm“M 





P2 

7665.5745 

7665.908315 

0 

7665.285014 

5 

PI 

7709.6839 

7710.017692 

0 

7709.3928 

0 

RO 

7788.8562 

7789.190809 

0 

7788.5673 

0 

R1 

7823.8212 

7824.162003 

5 

7823.525966 

5 

R2 

7855.6428 

7855.983584 

5 

7855.343839 

2 

R3 

7884.2790 

7884.621637 

0 




TABLE 4 — Continued 


(c) 2-0 band 


Line 

Standard Position (cm 

Venus A 


Venus 5 



(±10~^ cm~M 

Central Position 

(cm“M 

Quality 

Factor 

Central Position 

Quality 

Factor 

P6 

5524.9996 

5525.235831 

4 

5524.7990 

0 

P4 

5577.3313 

5577.571096 

2 

5577.118853 

5 

P3 

5601.7659 

5602.009737 

0 

5601.552334 

2 

P2 

5625.0282 

5625.268599 

5 

5624.812772 

2 

RO 

5687-6504 

5687.894402 

4 

5687.433641 

3 

Rl 

5706.0940 

5706.338637 

5 

5705.877340 

3 

R2 

5723.3020 

5723.548570 

1 

5723.085991 

2 

R3 

5739.2627 

5739.508974 

5 

5739.043513 

5 

R4 

5753.9647 

5754.210921 

4 

5753.744724 

5 

R5 

5767.3970 

5767.643448 

4 

5767.176221 

5 

R6 

5779.5488 

5779.795288 

5 

5779.327061 

5 

R7 

5790.4098 

5790.656031 

4 

5790.187529 

2 

R8 

5799.9700 

5800.216745 

4 

5799.746375 

1 

R9 

5808.2197 

_ 

- 

5807 . 9968 

0 


TABLE 4— Continued 




(d) 2-0 band 




Line 

Standard Position (cm~^> 

Venus A 


Venus 5 



(±10~3 cm“^) 

Central Position 

Quality 

Factor 

Central Position 
<cm 

Quality 

Factor 

P5 

5547.8656 

- 

- 

5547.6511 

0 

P4 

5573,4171 

5573.655490 

2 

- 

— 

P3 

5597.8123 

5598.053033 

2 

- 

— 

P2 

5621.0379 

5621.278192 

1 

- 

— 

PI 

5643.0807 

5643.322226 

2 

5642.866099 

1 

RO 

5683.5662 

5683 , 809110 

4 

5683.349254 

4 

R1 

5701.9840 

5702.222899 

0 

5701.759798 

0 

R2 

5719.1689 

5719.413631 

2 

5718.950232 

2 

R3 

5735.1094 

5735.354644 

4 

5734.890669 

5 

R4 

5749.7940 

5750.039806 

5 

5749.574644 

5 

R5 

5763.2118 

5763.457811 

3 

5762.991033 

4 

R6 

5775.3521 

5775.597533 

3 

5775.131670 

4 

R7 

5786.2046 

5786.448285 

2 

5785.9822 

0 

R8 

5795.7594 

5796.008412 

0 
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TABLE 5 

VENUS DOPPLER SHIFTS CALCULATED FROM EPHEMERIS DATA 


Doppler Shift. Factor (x 10 


Spectrum Date 

Planetary Motions 

Earth. Rotation 

Total 

143 

14 Oct 73 

+4.2986 

-0.0155 


+4.2831 

149 

15 Oct 73 

+4.3035 

+0.0090 


+4-3125 

156 

22 Oct 73 

+4.3359' 

-0.0268 


+4.3091 

161 

25 Oct 73 

+4.3488 

-0.0205 


+4.3283 

Venus A 

average 

+4.3217 

-0.0134 


+4.3083 

226 

13 Jun 74 

-3.9204 

+0.1105 


-3.8099 

228 

13 Jun 74 

-3.9013 

+0.1056 


-3.7957 

230 

14 Jun 74 

-3.8823 

+0-1119 


-3.7704 

231 

15 Jun 74 

-3.8627 

+0.1112 


-3.7515 

23'4 

16 Jun 74 

-3.8426 

+0-1103 


-3.7323 

Venus 5 

average 

-3.8819 

+0.1099 


-3.7720 
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TABLE 6 

BEFERENCE HEIGHTS OF CHORDS ACROSS LINES IN A-H(l) MODEL 


Line 



Height (km) 




Cl 

<=2 

C3 

C4 

Cs 

RO 

71.7 

71.4 

70.2 

67.4 

67.8’ 

Rl 

71.1 

70.2 

69.1 

66.0 

67.1 

R2 

f 

71.1 

69,9 

69.0 

65.9 

67.1 . 

R3 

69.2 

68.4 

67,6 

66.4 

65.8 

R4 

68.7 

68.6 

68.1 

67.2 

65.2 

R5 

68.9 

68.7 

68.2 

67.2 

64.2 

R6 

68,8 

68.6 

68.1 

67.4 

63.6 

R7 

68.5 

68.3 

67.9 

67.3 

63.9 

R8 . 

68.3 

68.0 

67.8 

67.1 

63.9 

Venus A 
Average 

69.5±1.3 

69.1+1.1 

68.410.8 

66.8+0.6 

65,411.6 

RO 

71.3 

70.6 

69.5 

66.3. 

67.3 

Rl 

70.0 

69.2 

67.7 

65.9 

66.2 

R2 

66.3 

66.0 

65,9 

65.5 

65.3 

R3 

69.4 

68.5 

> 

67.2 

65.9 

66.1 

R4 

69.4 

68.7 

67.8 

66.3 

65.9 

R5 

69,2 

68.9 

68.2 

66.9 

65.2 ' 

R6 

69.0 

68.7 

68.3 

67.4 

63.7 

R7 

68.7 

68.4 

68.1 

67.3 

63.7 

R8 

68.4 

68.1 

67.9 

67.2 

63.8 

Venus 5 
Average 

69.1+1.3 

68,5+1,2 

67.8+1.0 

66.5+0.7 

65. ’211. 3 

Overall 

Average 

69.3±1.3 

68.811.2 

68.1+0.9 

66.7+0.7 

65.311.4 
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TABLE 7 

-■ VENUS WIND VELOCITIES CALCULATED FROM C^^O^® 2-0 BAND 


Doppler Shift Factor Wind Velocity 

(x 10 (m sec 


Spectrum 

Number of 

CO band 

Ephemeris 

V 

U 


Lines 





143 

14 

+4.2696 

+4.2831 

+ 40.5 

+ 36.0 

161 

14 

+4.2795 

+4.3091 

+146.3 

+124.4 

Venus A 

14 

+4.2726 

+4.3083 

+107.0 

+ 93.2 

226 

16 

-3.8601 

-3.8099 

+150.5 

-169*. 1 

234 

16 

-3.8171 

-3.7323 

+254.2 

-293.9 

Venus 5 

16 

-3.8406 

-3.7720 

+205.7 

-234.3 
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Fig. 1. — Measured pressure-induced lineshifts for HF per atm Ar. 
(a) 1-0 band. Circles: Guelachvili and Smith C1978) ; Triangles; Jaffe 

et al . (1965); Crosses; Oksengorn (1963). (b) 2-0 band. Circles; 

Guelachvili and Smith (1978); Crosses: Wiggins et al. (1970. 
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Fig. 2. — Measured lineshifts per atm Ar in the HCil 2-0 band. (a) H Cfi.. 

Circles; Guelachvili and Smith (1978); Triangles; Rank et al . (1960); Inverted triangles; 
Rank et al. (1964); Crosses: Hirshfeld et al. (I960). (b) H CJl. Circles; Guelachvili 

and Smith (1978); Crosses; Hirshfeld ^ al . (1960). 
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’ Fig. 4. — Self-induced lineshifts in the HCJl 2-0 
measured at low pressure. Circles: Triangles: H CJl 

(Guelachvili and Smith, 1978). Dashed line: Kimel et al . (1959). 
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Fig. 5. — Self-induced lineshifts in the HF 2-0 band 
measured at low pressure. Circles: 'Observed 'shifts 

(Guelachvili and Smith, 1978) . Dashed lines: Shifts cal- 

culated by Boulet et al. (1976) . 
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Pig. 6. — Ground-state rotational energy levels and 
near-resonant transitions for and HF (after Boulet 

et al. 1976) . 
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Fig. 7. — Lineshifts in the HCA 2-0 band perturbed by 20 
torr -HF. Circles: observed shifts; triangles; ob 
served shifts (Guelachvili and Smith, 1978); inverted -triangles 
calculated shifts (Boulet, 1977) • 
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Pig. 8. — Lineshifts in the HF 2-0 band perturbed by 10 
torr HCJl (Guelachvili and Smith, 1978) . 
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Fig. 9. — Measured lineshifts per atm N2 in the HCit 2-0 band. Circles; Guelachvili 
and Smith (1978); triangles: Rank ^ (1960a); Crosses; Shifts measured by Kirael et 

(1959) for 1 atm dry air. (a) (b) H^^CA. 








Fig- 10 — Measured lineshifts for HF perturbed by CO 2 . (a) 1-0 band. Circles 
Guelachvili and Smith (1978); Triangles; Shaw and Lovell (1969). (b) 2-0 band. 

(Guelachivili and Smith, 1978 ) . * 
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Fig. 11. — Profile of temperatures in the Venus atmosphere. 
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Fig- 12. — Measurements of Venus zonal winds. 
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Flg._16. Residual lineshifts for HF in the Venus spectra 
^ cm ^). (a) 1-0 band, (b) 2-0 band. 
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Fig. 18. — Variation of translational Doppler shift factor 
from 1971 to 1980 (Niehaus and Petrie, 1961). 
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Fig. 20. — Variation of brightness along parallels on the planetary 
disk. Solid curves: Lambert's Law reflection. Dashed curves: isotropic 
scattering . 
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Fig. 22. — Pressure, temperature and velocity profiles for 
Venus model atmospheres. 
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Fig. 23 » — Absorption coefficient for significant 
for the'.H^^C^. R4 Line. • (a) Venus A 


levels in the median model model atmosphere 
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(b) Venus 5 
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Fig. 25. — Calculated asymmetric profile for the R4 line in the median 

model atmosphere for Venus 5. 
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Fig. 28. — Calculated lineshifts in the R-branch for line profile “D” 

in the A-H(2) model atmosphere. (a) Venus A. (b) , Venus 5. . . * 
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Pig. 30. — Weighted-average lineshifts in the HCi!- 2-0 band for the averaged 
planetary spectra. (a) Venus A. (b) Venus 5. 
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DRIGINAr PAGE IS 

Pig. 31. — Venus zonal velocity profilli^e^^S^k^S^St^ 
interferometric spectra. Triangles: Venus A; Circles: Venus 

143 (first day); Rhombi: Venus 161 (last day) . Diamonds: 

Venus 5; Squares: Venus 226 (first day); Crosses: Venus 234 
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